ECO 341 

Mathematics for Economists 


i 




11 



Ibadan Distance Learning Centre Series 


ECO 341 

Mathematics for Economists 


By 

Kareem Olayinka Idowu 

Department of Economics 
University of Ibadan 

and 

Olubiyi Ebenezer Adesoji 

Department of Economics 
University of Ibadan 



Published by 
Distance Learning Centre 
University of Ibadan 



© Distance Learning Centre 
University of Ibadan, 
Ibadan. 


Al! rights reserved. No parí of this publication may be reproduced, stored 
in a retrieval system, or transmitted in any form or by any means, 
electronic, mechanical, photocopying, recording, or otherwise, without the 
prior permission of the copyright owner. 

First Published 2009 


ISBN 978-021-402-X 


General Editor: Prof F. O. Egbokhare 
Series Editor: Olubunmi. I. Adeyemo 


Typesetted @ Distance Learning Centre, University of Ibadan 



Table of Contents 


Page 


Vice-Chancellor's Message. ... ... ... ... vi 

Foreword. ... ... ... ... ... ... vii 

General Introduction to the Course ... ... ... ... viii 

Lecture One: Linear Models and Matrix Algebra. 1 

Lecture Two: Matrix Algebra 2(Transposes Determinants 

and Inverses) ... ... ... ... 7 

Lecture Three: Multivariate Differential Calculus... 7S. 26 

Lecture Four: Optimization Problems... ... ... 39 

Lecture Five: Function Approximation. ... 60 

Lecture Six: Integral Calculus (Indefinite Integral) ... 39 

Lecture Seven: Differential and Difference Equations ... 79 

Lecture Eight: Optimal Control Theory... ... ... 92 


v 



Vice-Chancellor’s Message 


I congratúlate you on being part of the historie evolution of our Centre for 
External Studies into a Distance Learning Centre. The reinvigorated 
Centre, is building on a solid tradition of nearly twenty years of Service to 
the Nigerian community in providing higher education to those who had 
hitherto been unable to benefit from it. 

Distance Learning requires an environment in which learners 
themselves actively participate in constructing their own knowledge. They 
need to be able to access and interpret existing knowledge and in the 
process, become autonomous learners. 

Consequently, our major goal is to provide full multi media mode of 
teaching/learning in which you will use not only print but also video, 
audio and electronic learning materials. 

To this end, we have run two intensive workshops to produce a fresh 
batch of course materials in order to increase substantially the number of 
texts available to you. The authors made great efforts to inelude the latest 
information, knowledge and skills in the different disciplines and ensure 
that the materials are user-fríendly. It is our hope that you will put them to 
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Foreword 


The University of Ibadan Distance Learning Programme has a visión of 
providing lifelong education for Nigerian citizens who for a variety of reasons 
have opted for the Distance Learning mode. In this way, it aims at 
democratizing education by ensuring access and equity. 

The U.I. experience in Distance Learning dates back to 1988 when the 
Centre for Extemal Studies was established to cater mainly for upgrading the 
knowledge and skills of NCE teachers to a Bachelors degree in Education. 
Since then, it has gathered considerable experience in preparing and 
producing course materials for its programmes. The recent expansión of the 
programme to cover Agriculture and the need to review the existing materials 
have necessitated an accelerated process of course materials production. To 
this end, one major workshop was held in December 2006 which have 
resulted in a substantial increase in the number of course materials. The 
writing of the courses by a team of experts and rigorous peer review have 
ensured the maintenance of the University’s high standards. The approach is 
not only to emphasize cognitive knowledge but also skills and humane valúes 
which are at the core of education, even in an ICT age. 

The materials have had the input of experienced editors and illustrators 
who have ensured that they are aecurate, current and learner friendly. They 
are specially written with distance leamers in mind, since such people can 
often feel isolated from the community of learners. Adequate supplementary 
reading materials as well as other information sources are suggested in the 
course materials. 

The Distance Learning Centre also envisages that regular students of 
tertiary institutions in Nigeria who are faced with a dearth of high quality 
textbooks will find these books very useful. We are therefore delighted to 
present these new titles to both our Distance Learning students and the 
University’s regular students. We are confident that the books will be an 
invaluable resource to them. 

We would like to thank all our authors, reviewers and production staff 
for the high quality of work. 


Best wishes. 



Professor Francis O. Egbokhare 

Director 
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General Introduction 


This course is an intermedíate mathematics for economists built on the 
basic mathematical foundation set out in the preceding classes, that is, 
ECO 104 and ECO 204. However, to facilítate reading, the first few 
lectures were deliberately prepared in form of revisión of what were 
taught in the earlier classes. This notwithstanding, it is expected that 
students are equipped with basic mathematical understanding before 
proceeding to read this text material. 

Each lecture discusses basic mathematical principies as well as their 
applications to economics. Also, after the discussion in each lecture, 
students are expected to practice some simple exercises. Since the essence 
of teaching this course is to be able to handle economic principie, models 
and theories mathematically, most exercises were drawn from economics. 

It is our view that students will find this material useful particularly in 
other economic courses both basics and applied. 



LECTURE ONE 


Linear Models and Matrix Algebra 


Introduction 

From your knowledge of market equilibrium analysis, you will notice that 
for a two-commodity market, it is relatively easy to establish and solve for 
equilibrium of variables of interest. However, in the real world, markets 
are interwoven. For instance, the demands for beans may not only depend 
on its own price and price of palm-oil, but also on the prices of rice, 
pepper, and even yam. In particular, it is possible to have n-commodity 
market in which case, it will certainly be difficult to solve for the variables 
of interest simultaneously. In this ehapter, we shall present a method that 
is capable of handling this problem and the related ones. This method is 
called matrix algebra. Matrix algebra pro vides compact way of writing an 
equation system, even an extremely large one. Also, matrix algebra 
allows us to test for the existence of a solution by evaluation of a 
determinant. Furthermore, it gives a method of finding that solution, if it 
exists. Matrix algebra becomes more useful when we encounter 
comparative and dynamic analyses and also in handling complex (or 
simple) optimization problems. 

Objectives 

At the end of this lecture, you should be able to: 

1. arrange equation systems in matrix format; 

2. find the determinant of a simultaneous equation; 

3. compute the inverse of matrix; and 

4. apply Cramer’s rule in solving simultaneous equation system. 
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CONTENT 

A matrix is a rectangular array of numbers, parameters, or variables. The 
members of the array, referred to as the elements of the matrix, are usually 
enclosed in brackets or sometimes in parenthesis or with double vertical 
lines. 

Consider a general system of m linear equations in n variables (xi, X 2 , 
X 3 ..,.x n ). Such linear equation can be written as follows: 

a n xi + ai 2 x 2 + ai 3 x 3 +.+ ai n x n = Vi 

a 2 1 x 1 + a 22 x 2 + a 23 x 3 +.+ a 2 n x n = v 2 

a 3 ixi + a 32 x 2 + a 33 x 3 +.+ a 3 n x n = v 3 . 1.1 


a m iXi + a m 2 x 2 + a m 3 x 3 +.+ a mn x n — v m 

In equation (1.1), the variable xi appears only within the leftmost column 
or in general, the variable Xj (j = 1, 2, 3... m) appears only in the jth 
column on the left side of the equal sign. The double-scripted parameter 
symbol a¡j represents the coefficient appearing in the ith equation and 
attached to the jth variable. As an example, a 32 is the coefficient in the 
third equation, attached to the variable x 2 . The parameter v¡ represents the 
constant term in the ith equation. Therefore, equation (1.1) contains three 
ingredients in equation system. These are: 

• a set of coefficient, a¡j 

• a set of variables, xi, x 2 , x 3 ,., x n 

• a set of constant terms, vi, v 2 , v 3 ,. , v m . 

it is possible to arrange this as arrays and label them respectively as 
matrices A, X and d as in: 
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j $12 $13 

$ 2 i $22 $23 

$31 $32 $33 


a ml a m2 a m3 

Let us give an example to illustrate further. Assuming we have a linear- 
equation system as follows: 

3xi + X 2 + 5x3 = 6 
5xi - 3x2 + X 3 = 15 
7xi - X 2 +3X3 =11 

We can use the idea of equation 1.2 to write the above equation in matrix 
form as follows 



What you should note is that the coefficient matrix both in equations 1.2 
and 1.3 are separated not by commas but by blank spaces only. This is 
how we represent linear-equation system in matrix form. 

Commutative, Associative and Distributive Laws of Matrix 

In ordinary scalar algebra, the additive and multiplicative laws obey the 
commutative, associative and distributive laws as follows: 

Given matrices A, B and C, 
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• commutative law of addition: a + b = b + a 

• commutative law of multiplicadon: ab = ba 

• associative law of addition: (a + b) + c = a + (b + c) = b + (a + c) 

• associative law of multiplication: (ab)c = a(bc) 

• distributive law: a(b + c) = ab + ac 

However, in matrix operations, most, but not all of these laws are 
applicable. The commutative law and distributive laws of addition hold in 
all cases. For instance, given matrices A, B and C, the commutativé law of 
addition says that 

A + B = B + A 

For associate law, (A + B) + C = A + (B + C) 

Thus (A + B) + C = (B + C) + A 

Matrix Multiplication 

Matrix multiplication is not in all cases commutative, for instance, given 
matrices A, B and C, AB ^ BA. This holds in all cases except when A is a 
squared matrix and B is an identity matrix, that is 

AB = BA if A is a squared matrix and B = I 
Where I is identity matrix. Another case is when matrix A is the inverse 
of matrix B, that is 

AB = BA when A = B' 1 

Another case where commutative law applies is in the case of scalar 
multiplication. That is given scalar v, and matrix A, then 
Av = vA. 

Associative law of multiplication 

For matrices ABC to be conformable for multiplication, it must be the 
case that each adjacent pair of matrices is satisfied. That is, if A is of mXn 
dimensión and C is vXk dimensión, then conformability requires that B is 
nXv 
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ABC 
mXn nXv vXk 


If the conformability condition is met, the associative law States that any 
adj acent pair of matrices may be multiplied out first, provided that the 
product is duly inserted in the exact place of the original pair. 

Distributive law of Matrix Multiplication 

Matrix multiplication is also distributive. In the case of distributive law, 
the conformability conditions for addition as well as for multiplication 
must be observed, i.e., given matrices A, B and C, 

A (B + C) = AB + AC (premultiplication by A) 

(B + C) = BA + CA (post multiplication by A) 

Identity Matrices and Nuil Matrices 

Identity matrix is defined as a square matrix with l’s in its principal 
diagonal and O’s everywhere. Principal diagonal is the line running from 
North-West (extreme up left) to South-East (extreme down right). We 
denote identity matrix by I. If we want to indicate identity matrix, we may 
write I2, I3 meaning that the identity is of two dimensión and three 
dimensión respectively. In general, I n (n = 1, 2, 3... k) means identity of 
dimensión n. 

One important role played by identity matrix is that it represents 
number 1 in scalar algebra. As we know, for any number a, we have l(a) 
or a(l) = a. Also for any matrix A, we have this relation: 

IA = AI = A.1.5) 

Identity matrix has a special feature in matrix multiplication. For instance, 
it is possible to insert or delete an identity without affecting the matrix 
product since AI = IA = A 
Now consider the associative law of the form 
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AI B (AI)B = AB 

mXn nXn nXk = n mXn nXk 


which shows that the presence or absence of I does not affect product 
operation? 

An identity matrix remains unchanged when it is multiplied by itself 
any number of times. Any identity matrix with such property is called an 

idempotent matrix. 

Idempotent matrix = AA =A 
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LECTURE TWO 


Matrix Algebra 2(Transposes Determinants 
and Inverses) 


Introduction 

Transpose of a matrix is obtained by interchanging the row to become the 
column. For instance, if we have matrix A and B such that: 



00 

m 

-9' 


'A 6" 

A = 



B = 



^6 0 

7 J 


vi 8 y 


We can interchange the row and columns to get 


( 3 

63 


(4 I a 




00 

0 

ii 

«5 


1-9 

7 , 


v 6 8 y 


By definition, if A is mXn dimensión, then its transpose, A, will be nXm 
dimensión. 

Consider the following matrix: 


A = 


"1 0 
0 3 

, 4 7 


7 


7 





Its transpose is: 


1 0 4 
Á= 0 3 7 
v 4 7 2 

which is the same as the original matrix. This class of squared matrix is 
known as symmetric matrix. Another example of this type of matrix is 
the identity matrix I, because its transpose, I, is actually I. Thus, we must 
bear it in mind that all identity matrix are symmetric matrix, while all 
symmetric matrix are not identity matrix. 

Determinants and Nonsingularity 

Determinant is a single number or a scalar. Determinant is applicable only 
to squared matrix. Let us consider a 2X2 matrix as in 



The determinant: I Al is found by taking the difference of the product of 
the two elements of the principal diagonal (ad) and the product of the two 
elements of the principal diagonal (cb). Thus 

I Al = ad-bc.2.4 


Example 



= 8-15 = -7 


(2.5) 
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Concept of nonsingularity 

Now consider the following 2X2 matrix: 


A= 


"4 

4^| 


r -3 3 N 



or\ 


v4 

4; 


v 6 -<C 


( 2 . 6 ) 


The determinant of each of the matrices is 0. This type of matrix is said to 
be singular simply because its determinant is zero. A matrix is said to be 
nonsingular when its determinant is not equal to zero (it could be negative 
or positive but not zero). The feature of a non-singular matrix is that there 
exists linear dependence between at least two rows or two columns. In 
short, if I Al = 0, it is a singular matrix and there exists linear dependence 
among equations. The implication of singular matrix is that there will not 
be a unique solution. Conversely if IAI ^ 0, it is a non-singular matrix and 
there is no linear dependence among the equations. The implication of 
this is that a unique solution exists. 

The rank of a matrix allows for a simple test of linear dependence and 
henee singular matrix. The rank p of a matrix is defined as the máximum 
number of linearly independent rows or columns in the matrix. Assuming 
a square matrix of order n, if p (A) = n, A is nonsingular and there is no 
linear dependence. If p (A) < n, A is singular and there is linear 
dependence. In 2.5, the matrix has a full rank because the rank of A is 2, 
written as p(A) = 2, and therefore it is easy to see that it is nonsingular. 
However, in 2.6, we see that row 2 and column 2 are twice times row 1 
and column 1 respectively. Henee p(B) = 1 and since p(B) < n (where n = 
2), then B is singular and linear dependence exists between rows and 
column. In this case there is no unique solution to the equation system. 

We now extend our discussion on determinant to 3X3 dimensión. 
Consider matrix A with the following elements: 


^ a b c' 
d e f 
K g h i , 


(2.7) 
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There are several ways of getting the determinant of this matrix. For our 
own purpose, we shall treat two. 

Steps to get determinant using deleting method 

a. assign + and - to the first row starting with + 

b. take the first element of the first row (a) and mentally delete the 
row and column in which it appears. 

c. Put - and take the second element in the first row (b) and delete its 
corre sponding row and column 

d. Put + and take the third element in the first row (c) and delete its 
related row and column. 

If this is done very well, and you apply the definition of determinant 
stated earlier, you will get something like this: 

a(ei-hf) - b(di-gf) + c(dh-ge)...( 2 . 8 ) 

Example: find the determinant of this matrix: 

Í1 2 6' 

A= 5 4 3 .(2.9) 

v 0 9 8 J ^ 

IAI =l[4(8)-9(3)]-2[5(8)-0(3)]+6[5(9)-0(4)l = 195 

The second method is by adding the producís of the principal diagonal, its 
minors and subtracts it from the product of the off diagonal elements and 
its minors from it. Let us apply this method to equation 2.9. 

"1 2 ó') 

A= 5 4 3 

[o 9 8 J 

IAI = [(Ix4x8)+(5x9x6)+(0x2x3)-[(0x43x6)-(9x3xl)-(8x5x2)l 
= 32+270+0-0-27-80 
= 302-107 

= 195 
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We see that we get the same answer. Note that the second method is 
applicable only to 3X3 matrix. If you have a matrix that has more than 
3X3 dimensión, the second method is irrelevant. 

Properties of a Determinant 

There are seven properties of a determinant. The se are: 

i. Adding or subtracting any non-zero múltiple of one row (or 
column) from another row (or column) will have no effect on the 
determinant. 

ii. Interchanging any two rows (columns) of a matrix will change the 
sign, but not the absolute valué of the determinant 

iii. Multiplying the elements of any row or column by a constant will 
cause the determinant to be multiplied by the constant 

iv. The determinant of a triangular matrix, i.e. a matrix with zero 
elements everywhere above or below the principal diagonal, is 
equal to the product of the elements on the principal diagonal 

v. The determinant of a matrix equals the determinant of its transpose 

vi. If all the elements of any row(column) are zero, the determinant is 
zero 

vii. If two rows (columns) are identical or proportional, i.e. linearly 
dependent, the determinant is zero. 

Matrix Inversión 

An inverse matrix A" 1 is a unique matrix satisfying the relationship: 

AA' 1 = 1 = A 1 A.2.10 

That is, if you multiply matrix A by its inverse, you will get an identity 
matrix. It then follows that in verse matrix in linear algebra performs the 
same function as the reciprocal in ordinary algebra. The formula for 
deriving the inverse is given by: 

A" 1 = -^adjunctA 
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Adjunct A is found by the transpose of cofactor of matrix A. The cofactor 
of matrix A is given by this formula: 


C¡| = (-1)' 1 j: |Mí|I .2.12 

Where Cy is the cofactor of minor My. My is the determinants of minors 
when row i’ is and column j is deleted. 

For illustration, let us begin with finding the inverse of a 2x2 matrix. 
Consider the matrix given in 2.3, in which the determinant is derived by 
equation 2.4. By equation 2.11, we have found IAl. What is left is to find 
the adjunct of A. To do this, we first find the determinants of the minors. 
These determinants are: 

Mn = d M 12 = c M 21 = b M2í^a 
Henee we have 



Having got the adjunct, we can now invoke equation 2.11 to get our 
inverse as follows 



ad -be 
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It is easy to demónstrate that AA" 1 = I. To do this, premultiply the inverse 
by the original matrix A as in: 


' d 

a b ad-bc 
v c d J —c 
V ad-bc 


-b 

ad-bc 

a 

ad-bc 


/ ad cb -ab ab ' 

-- 1 - 

ad-bc ad-bc ad-bc ad-bc 

cd cd -be ad 

-- 1 - 

V ad-bc ad-bc ad-bc ad—be) 


1 0 
0 1 


Solving Linear Equations using matrix Inversión 

An inverse matrix can be used to solve matrix equations. If A n x n X n xi B n xi 
and the inverse A" 1 exists, multiplication of both sides of the equation by 
A' 1 , and following the laws of conformability gives: 

A n.XnA n X]XnX | — A nX n B n XI 

From equation 2.10, A" A = I thus 
InXnX n xi — A nXnBnXl 


Also consider our equation 1.4 the matrices can be written as 


A = xv 

Since we are looking for valúes of x’s then we can write this as 

- = x = A _1 v.(2.13) 

d 

Thus, we have a powerful tool for solving system of equation. For 
instance, the system of equation presented in equation 1.4 can be solved 
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simultaneously using equation 2.10. for convenience, let us rewrite 
equation 1.4 derived from equation 3: 



A x = v 

Ax = v .-. x = A v.(2.14) 

The determinant is: 

lAl = 3(-9-l) -1(12+7) +5(-4+21) 

-30-19+85 

36 ^ 0.(2.15) 

Next we find the determinants of the minors: 



IMnl = -10 IMi 2 I = 19 IM13I = 17 

\Ú 2 i\ = 8 IM 22 I = -26 IM 23 I = -10 

IM31I = 14 IM32I = -23 IM33I = -13 
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Thus M¡j = 


-10 19 17 " 

8 -26 -10 
14 -23 -13, 


The cofactor then becomes: 


C = 


'-10 -19 17 " 

-8 -26 10 
v 14 23 -13, 


C' = transpose of C = adjunct of C 




'-10 

-8 

14 

Adj A 

= 

-9 

-26 

23 



V 17 

10 

-13 

The inverse of matrix A 

A" =4. 

'-10 

-8 

14 

-9 

-26 

23 

36 





l 

17 

V 

10 

-13 


- 

10 

-8 

14 “ 


36 

36; 

36 



-9 

-26 

23 


36 

36 

36 


17 

10 

-13 


36 

36 

36 
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X 


-10 

-8 

14 

36 

36 

36 

-9 

-26 

23 

36 

36 

36 

17 

10 

-13 

36 

36 

36 


6 

15 

11 


after simplifying, we find out that 


Xj = -0.72, X 2 = -7.56 X3= 3.02 

Solving Linear Equation using Cramer’s rule 

Another approach used to find solution to linear equation system is the use 
of Cramer’s rule. Cramer’s rule provides a simplified method of solving a 
system of linear equation through the use of determinants. Cramer’s rule 
is given by equation 2.13. 


I Al 

where a¡ is the ith unknown variable in a series of equations, IAI is the 
determinant of the coefficient matrix and IA¡I is the determinant of a 
special matrix from the original coefficient matrix by replacing the column 
of coefficients of x¡ with the column vector of constants. We can apply 
Cramer’s rule to solve equation 4. 

The determinant of A has been found in equation 2.2. Thus, we need to 
solve for IAiI, IA 2 I, and IA 3 I as follows: 


Ai = 


6 

15 

11 


1 5 

-3 -1 
-1 3 
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IAil = 6(-9-l) -1(45+11) +5(-15+33) 

-60-56+90 

-26 

'3 6 5 ' 

A 2 = 4 15 -1 

7 113 

|A 2 I = 3(45+11) -6(12+7) +5(44-105) 

3(56) -6(19) +5(-61) 

-251 

*3 1 6" 

A 3 = 4 -3 15 

7 -1 11 

|A 3 | = 3C-33+15) -1(44-105) +6(-4+21) 

-54+61+102 = 163-54 
109. 

I leave getting ai, a 2 , and a 3 as an exercise (Hint: note that á¡= \A¡\AA\, 
ü 2 = IA 2 I/IAI; Ü 3 — \A 3 \AA\) 

Application of the use of Matrix inversión and Cramer’s rule to 
Economics 

A case of a 3-commodity Market 

At the beginning of chapter one, we pointed out that the demand for beans 
may not only depend on its own price, but also depends on the prices of 
rice and even kerosene. Let us assume that the demand for and supply of 
quantities of beans, rice and kerosene is given as follows: 
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Qdb = 4 - 3Pb + P r + 2P k Q S b = -3 +2Pb 

Qdr = 2 + Pb~ P r + Pk Qsr = -2 + P r 

Qdk = 1 + 2Pb + P r - 2Pk Q sk = -1 + 2P k 

Where Qj, and Q s ¡ (i = b, r, k) are quantities demanded and supplied of 
beans, ríce and kerosene, while P¡ (i = b, r, k) are prices of beans, rice and 
kerosene respectively. Assuming that we are interested in finding the 
equilibrium prices of quantities of the three markets, we can use matrix 
inversión or Cramer’s rule to do this. 


The first step to take is to establish the equilibrium condition and 
simplify the linear equation such that it becomes a system of equation. 
The equilibrium condition is given by: 

Qdb - Qsb = 0 .. < v / ..2.17 

Qdr - Qsr =0 .; .. 2.18 

Qdk-Qsk=0 . y. .2.19 

Thus 4 - 3P b + Pr + 2Pk - [-3 +2P b ] = 0 

4 - 3P b + Pr + 2P k + 3 - 2P b = 0 

7 -5 P b +P r + 2P k =0 

-5 P b + P r + 2P k = -7. 

.2.20 

2+ Pb-Pr+.Pk- [-2 + P r ] = 0 

2 + Pb - P r + Pk + 2 - P r = 0 

4 + Pb - 2P r + Pk = 0 

P b 7 2P r + P k = -4. 

.2.21 

1 + 2P b + P r - 2P k -[-l + 2P k l 

1 + 2P b + P r - 2P k + 1 - 2P k = 0 

2 + 2P b + Pr -4 Pk = 0 

2P b + P r -4 P k = - 2. 

.2.22 


Equations 2.20, 2.21 and 2.22 form the system of linear equation for 
which equilibrium conditions are satisfied. As can be seen, we have 3 
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equations in 3 unknowns. To solve for the unknowns using matrix 
inversión, we first write equations 2.20, 2.21 and 2.22 in matrix format as 
follows 


'-5 1 2 ^ 


ÍPb' 


(-T 

1 -2 1 


Pr 

= 

-4 

,2 1 -4) 


KPk, 


l-2 y 


We can then proceed by first finding the determinant of the matrix of 
coefficient, the cofactor, the adjunct and the inverse. 

The determinant of the matrix of coefficient, cali it A is calculated thus: 

I Al = -5(8-1)-l(-4-2)+2(1+4) 

= -5(7)-l(-6) +2(5) 

= -35+6+10 
= -9 

Then we use equation 2.10 to solve for the equilibrium prices. Having got 
these prices, we then substitute the equilibrium valúes for the prices in 
demand and supply equation to get the equilibrium quantities. 

We can also use Cramer’s rule to solve for equilibrium prices by evoking 
equation 2.13. Do this as an exercise. 

A Case of National Income 

Consider the national income equation of the form: 


Y = C + I.2.23 

C = a + (3Y.2.24 

I = lo.2.25 


Where Y, C and I are national income (which equals aggregate planned 
expenditure), planned prívate consumption expenditure and planned 
prívate investment respectively. Again, a and (3 are autonomous 
consumption and marginal propensity to spend respectively. As can be 
seen, there are two endogenous variables in the model which we should 
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solve for. These endogenous variables are Y and C. We can use matrix 
inversión or Cramer’s rule to solve for the equilibrium Y (Y*) and C (C*). 

To do this, substitute 2.24 into 2.23 and re-arrange equations 2.23 and 
2.24 so that the endogenous variables together with their coefficients are 
on one side while the rest are on the other side as follows: 


Y - C = I 0 
- (3Y + C = a 

We can now write this in matrix form as follows: 

We thus have matrix of coefficient A, matrix of variables X and matrix of 
constants V. It is straightforward to verify that the determinant of matrix 
A, IAl, is (l-(3). The cofactor and the adjunct, respectively, are 


C = 


ri 

vi 


A" 

ij 


AdjA = 


íl 


n 

K 


A' 1 = 


1 


1 -b 


1 1 
A 1 


Since X = 
Then 
X = 


Ai 


\Cj 


A _1 V 


1 

“1 f 

V 

1 -b 

_A i. 

a 
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.2.26 


Henee, Y* 


1 


i-A 


(lo + a) 


C*= — 1 — (a+pl).2.27 

1 -p 

We can verify that this result is the same if we use the traditional 
substitution method. 

A Case of IS-LM relation in a closed economy 

Now, consider an economy with two specific sectors - the real sector and 
the monetary sector. The product market (which constitutes the real 
sector) has the following equilibrium condition and behavioural equations 

Y = C + I + G (equilibrium condition in the product market).2.28 

C = a + pY d (private planned consumption expenditure behavioural equation)..2.29 

Y d = Y -T (disposable income behavioural equation).2.30 

T = tY (tax revenue behavioural equation).2.31 

I = d - ei (private investment expenditure behavioral equation) .. .2.32 
G =Go (government final consumption expenditure behavioural 
equation) ..v.í... 2.33 

The endogeneous variables are Y, C, I and i (the real interest rate). 
The endogeneous variable is Go while a, d, e, [), and t is structural 
parameters. 

In the second sector, the monetary sector, the equilibrium condition 
and behavioural equations are specified below: 

Md = M s (equilibrium condition).(2.34) 

Md = kY - L¡ (money demand behavioural equation).(2.35) 

M s = Mo (money supply behavioural equation).(2.37) 


We see that money supply is exogeneous; henee, it is the high- 
powered money. The variables k and L are structural parameters of the 
money market. Using the equilibrium condition stated in equation 2.34, 
equations 2.35 and 2.36 become 
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M 0 = kY - L¡ 


2.38 


Equation 2.33 gives the equilibrium in the money market or the LM 

equation. 

In the product market, 

C = a + (3(l-t)Y.2.39 

Using equations 2.29, 2.39, 2.32 and 2.38, the two sectors produce the 
foliowing system of equations: 

Y -C-I = G 0 

(3( 1 -t)Y - C = - a 

I + ei = d 

kY -Li = M 0 

we have 4 equations in 4 unknowns. We can further reduce the system of 
equation to 2 equations in 2 unknowns. But to demónstrate the ability of 
matrix in handling equation of this nature, we deliberately leave it at 4 
equations. It is time to arrange the equations in matrix form and this is 
done below. 

1 - 1-1 0 Vf] ÍG 0 ' 

Z?(1-0'K-1 0 0 C -a 

H = .2.40 

0 0 1 el d 

K 0 0 -L)\i) v M 0y 

We can then proceed to finding the inverse of matrix of coefficient and so 
on. It is also possible to use Cramer’s rule to determine the equilibrium 
income, consumption, investment and interest rate. 

A case of Production Activities (The Leontief Input-Output Model) 

Assuming that we have three sectors in an economy producing goods and 
Services to be consumed by the population, let us assume that the 
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production technology in each sector is such that it uses some of the 
output produced in other firm (together with its own output) as input. This 
implies that the production of the three sectors depends on each other. 
The correct output demands are those which satisfy technical input-output 
relationships. The question that we ask is what is the output demand of 
each sector? This is the issue tackled by Leontief input-output model. 

The problem posed by input-output analysis can be solved with the 
aid of system of simultaneous equations, which can be conveniently 
handled by matrix algebra. If the first sector is to produce an output just 
sufficient to meet the input requirements of the three sectors as well as the 
final demand (final consumption), its output level, say Xi must be of the 
form: 

Xi = aiiXi+a^Xi+aoXa+di 

Where di denotes the final demand for its output and a¡jXj represents the 
input demand from the jth industry. In the same vein, the output levels of 
the other two sectors are 

X2 = a 2 iXi+a22X2+a23X 3 +d2 
X 3 = a3iXi+a32X2+a33X3+d32 

Which can be written in matrix form as? 


po 


^ Cl 11 $12 $13 




r d x ' 

x 2 

= 

$ 2 i $ 2 2 ¿*23 



+ 

d 2 



^£¿31 $32 $33 j 




1^3 y 


This can be written in compact forms as 

X = aX+d.2.41 

Where X is a vector of output, a is a vector of technical coefficients 
and d is a vector of final demand (external demand). To find the level of 
total output (intermedíate) needed to satisfy final demand, we can solve 
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for X in terms of the matrix of technical coefficients and the column 
vector of final demand as follows 


X-aX = d 


2.42 


X(l-a) = d 


X(I-a) = d.2.43 

X = (I-a) _1 d.2.44 


Where I is the identity matrix, we can write Equation 2.44 as follows 


x 2 


( í 


W 


1 0 0 a 
0 1 o 
0 0 1 


\ a 3 i 




a 32 a 33 J) 


d 2 

\d-3j 




— a n 

— a l2 

a \3 

x 2 

= 

— ü 2 \ 

1 — ü 22 

—«23 



l ~ a 3 i 

~ a 32 

1 «-33 j 


-i 


d\ 


.(2.46) 

Equation 2.46 is called input-output formula for a three-sector model 


\®3j 


.2.45 


Example 

Consider an economy consisting of three sectors, namely agriculture (A), 
manufacturing (F) and monetary sector (M). The production of agriculture 
uses 0.3 of its own output, 0.4 share of manufacturing and 0.1 share of 
monetary sector. The production in manufacturing takes place with 0.5 of 
its output share, 0.2 of agriculture’s share and 0.6 of monetary sector 
output. Finally, for monetary sector, 0.1 share of agriculture, 0.3 share of 
manufacturing and 0.1 of own outputs are needed. Given that the total 
final demand for each, 20, 10 and 30 respectively, determine the total 
demand for each sector. 


Set up the problem in system of equation as follows 
A = 0.3A+0.4F+0.1M = 20 
F = 0.5A+0.2F+0.6M = 10 
M = 0.1A+0.3F+0.1M = 30 


24 








Let A, F and M be Xi, X 2 and X 3 respectively. Thus applying equation 
2.46, we have 




( 0.7 

-0.4 

- 0 . 1 ) 

-1 

f20' 


= 

-0.5 

0.8 

- 0.6 


10 

l*3 j 


i-o-i 

-0.3 

0.9 J 


W 


First solve for the inverse of the technical coefficient matrix (1-a )" 1 
This yields the result: 


(I-a ) 4 


1 

0.151 


0.54 

0.39 

0.32 

0.51 

0.62 

0.47 

0.23 

0.25 

0.36 


By using equation 2.39 we find out that 


X = 


1 


0.151 


1 


0.151 


0.54 0.39 0.32 
0.51 0.62 0.47 
v 0.23 0.25 0.36 

^24.3 
30.5 
v 17.9 


20 

10 

v 30 y 


\ 

f 160.93' 

= 

201.99 

J 

^ 118.54 y 


Thus, the total demand for 
Agriculture sector = 160.93 or 161 
Manufacturing sector = 201.99 or 202 
Monetary sector = 118.54 or 119 
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LECTURE THREE 


Multivariate Differential Calculus 


Introduction 

Earlier discussion at the lower class on differential calculus was based on 
single variable. In comparative-static analysis, however, it is usual to 
come across a situation in which we have several variables and 
parameters, so that the equilibrium valué of each endogenous variable may 
be a function of more than one variable or parameter. To this end, we need 
to learn how to find the derivative of functions having more than one 
variable. This chapter deais with comparative-static models with more 
than one variable. 

Objectives 

At the end of this lecture, you should be able to: 

1. find the rate of change of any parameter in a multivariate model; 

2. differentiate between constrained and unconstrained optimization; 

3. identify and operationalize a set of rules that can be used to solve 
unconstrained and constrained optimization problems in a 
multivariate model; and 

4. analyze comparative static of any economic models. 

Pre-Test 

1. Given Y = f(a,b) = (a+2)(3a+7b), find the pardal derivatives of the 
function. 

2. Given the demand and supply equations: Qd = a - bP; Q s = -c - dP, 
find the equilibrium prices and quantities. Find the rate of change 
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in price with respect to parameters a, b, c and d. interpret your 
result. 

CONTENT 
Partial Derivative 

Consider the function: 

Y = f(zi, z 2 , z 3 ,., z n ).,3.1 

Where z¡ (i = 1, 2, 3... n) are all exogeneous variables such that each can 
vary by itself without affecting others. If z\ undergoes a change, Azi say, 
while others remain fixed, there will be a corresponding change in y, 
namely Ay. Thus, the process of examining a change in y with respect to a 
change in z\, holding other variables constant, is called partial 
differentiation. In particular, partial differentiation differs from other 
differentiations because we must hold (n-1) exogenous variable constant, 
while allowing one variable to vary. To differentiate partial derivative 
sign from other signs, it is conventional to use “curved d”, that is d. 

Illustration 1: consider the function y = f(zi, z 2 ) = 3z¡ z 2 

Let us assume we want to consider a change in y as a result of a change in 

z\. We will treat z 2 constant by mentally bracketing it with the coefficient: 

o 2 3 

y = 3 z¡z 2 

^ = yi = (3z 3 2 )^-(z¡) 
oz¡ azj 

= 3z 3 2 .2 = 6z! 

To examine the change in y with respect to change in z 2 , we treat z\ as a 
constant and follow the procedure above: 

dy 2, d 3 

— =y 2 =(3 Zl ).—(z 2 ) 

az 2 dz 2 

= (3zf).3z¡=9zh¡ 
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Illustration 2: find the partial derivative of the following function 
y = f(x,z) = 5x 3 - 3 x 2 z 2 + 7z 5 . 

dy d 3 2 n d 2 d 5 

— = —(5x )-(3z ). — (x ) + — (7z ) 

dx dx dx dx 

= 15x 2 - (3 z 2 ).2x + 0 (why 0?) 

— = 15x 2 - 6z 2 x 

dx 


dy 

dz 

dy 

dz 


= -¡~(5x 5 )-(3x 2 ).-^-(z 2 ) + 
dz dz 

= 0 - (3x 2 ).2z + 35z 4 

= -6x 2 z + 35z 4 . 


—(7z 5 ) 


Rules of Partial Differentiation 

1. Product rule: given y = f(xi, x 2 ).g(xi, x 2 ), 

j., . dg . . df 

= — + g(x lf x 2 ).^~ 


dx. 


dx t 


dx ; 



f(x v x 2 ). 



+ g{x v x 2 ). 



Example: given y = (3 xi+5)(2xi+6x 2 ) 

In this case, f(xi, x 2 ) = (3xi+5) 
While g(xi, x 2 ) = (2 xi+6x 2 ) 
Invoking product rule, 
dg_ _ 2 . _ 3 

3.Xj d.Tj 

=6 .df =Q 

dx 2 dx 2 

Henee, by equations 3.2 and 3.3, 


3.2 


3.3 
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P- = (3x 1 + 5).(6) + (2xi+6x 2 )(3) 
óx¡ 

= 6 xi +10 + 6 xi + 18 x 2 

= 12xi + 18x 2 + 10 

3xj 


= (3x 1+ 5)(6) + (2 x 1+ 6x 2 )(0) 

ox 2 

= 18xi + 30 + 0 

= 18xi+ 30 


2. Quotient Rule: y = ^ X| ’* 2 ; g(xi, x 2 ) ^ 0 

<?(*!, *2 

5y/5xi = g (*!, x 2 )-df / dx ] - \f(x ¡ ,x 2 ).3g /3x, 1 

C U<+ P +Jl : 

dy/dxi= g(x 1 ¿k 2 ).df/dx 2 - [f(x 1 ,x 2 ).dg/dx 2 ] 
[g(x 1 ? x 2 )] 2 


Example: let y = 


_ ( 6 x, + 7x 2 ) 


, find the partial derivatives of the function 


(5xj + 3x,) 

From the problem we know that: 
f(xi, x 2 ) = ( 6 xi + 7x 2 ) 
g(xi, x 2 ) = (5xi + 3x 2 ) 

thus, invoking quotient rule given by equations 3.4 and 3.5, 
— = 6 ; = 7 


3x, 


dx-, 


dg =5 . d S _ 3 

3xj 3 xt 


Henee, 

3v _ (5.x, +3 x 2 ).6-[(6.Xj + 7x 2 ).5] 
3xj (5Xj+3x 2 ) 2 
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_ 30.Xj +1 8x 2 - 30-Xj - 35x 2 
(5x, + 3 x 2 ) 2 

18x 2 -35x 2 
(5xj +3x 2 ) 2 

3y _ 17 x 2 

3x[ (5 xj + 3x 2 ) 2 

3y _ (5xj +3 x 2 ).7-[(6xj +7x 2 ).3] 

3x 2 (5xj+3x 2 ) 2 

_ 35xj + 2lx 2 -18xj - 2 lx 2 
(5xj +3 x 2 ) 2 

35Xj -18xj 

(5x t + 3 x 2 ) 2 
3y _ 17x t 

3x 2 (5xj + 3x 2 ) 2 r\'-' 

3. Generalized Function I Rule: Given y |f(x ¡, x 2 )] n 

= n[f(x h x 2)l # .3.6 

3x, ax, 

dy/dx 2 = n[f(x 1 ,x 2 )r 1 ^- .3.7 

ox 2 

Example: find the partial derivatives of the function y = (xi 3 + 7x 2 4 ) 5 
In the above example, 

<9f/<9xi = 3xi 2 ; df/dx 2 = 28x 2 3 
dy/dx! = 5[(xj 3 + 7x 2 4 ) 5 ‘l(3xr) 

3y/3xj= 15(3 xi 2 )(xi 3 + 7x 2 4 ) 4 

dy/dx 2 = 5[(xi 3 + 7x 2 4 ) 5 ‘l(28x 2 3 ) 

Sy/dXi = (140x 2 3 )( Xl 3 + 7 x 2 4 ) 4 
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Application to Economics 

1. Consider a simple one-commodity market model, in which case the 
demand and supply equations are given as follows: 

Q d = a - bP 
Q s = -c + dP 

the equilibrium Solutions are: 

p* CL + C 

b + d 

* _ ad -be 
Q “ b + d 

Now let us consider how a change in each of the parameters will affect the 
equilibrium price. To do this, we will partially differentiate P with 
respect to each parameter. Note that we have four parameters that 
determine the equilibrium price; therefore, we shall establish four partial 
derivatives. Which rule are we going to invoke? Quotient rule of course. 
Henee, we let f(xi, X 2 ) = a+c while g(x¡, X 2 ) = b+d. By equations 3.4 and 
3.5, 

dP* _ (b + d).l-[(a + c).Q] 
da (b + d ) 2 

b + d 
(b + d ) 2 


da b + d 

dP _ (b + d)0-[(a + c)ll 
~db~ (b + d) 2 

dP a + c 

~—= - 7 <0 

db (b + d)~ 

dP _ (Z? + £/)l-[(a + c)01 
3c (b + d) 2 
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= - >0 

b + d 

dP 8 _ (b + d)0-[(a + c)l] 
d d (b + d) 2 

dP _ a + c 
~dd~ (b + d) 2 < 

Henee cP/ca = cVIcc = 1/b+d > 0 

dvldb = tVIcd = - (a+c)/(b+d) 2 < 0 

That is, an increase in “a” (the intercept of demand) or an increase in 
“c” (the intercept of supply), will, holding other things constant, increase 
equilibrium price. Conversely, an increase in “b” (the slope of demand) or 
an increase in “d” (the slope of supply), will, holding other things 
constant, reduce equilibrium price. Since a and b appear in the demand 
equation, it means that changes in either of the two will affect the demand 
curve only. How will the effect take place? We know that a is 
autonomous while b is induced. Therefore, a change in equilibrium price 
with respect to a change in a will shift the demand curve parallel to itself, 
while a change in equilibrium price caused by a change in b will pivot the 
demand curve inward or outward depending on whether such change is an 
increase or a decrease in b. 

We also know that c and d are parameters appearing in the supply 
equation. Henee, a change in equilibrium price caused by changes in any 
of the two wiíl affect the supply curve only. Since c is an exogenous 
parameter, a change in equilibrium price with respect to c will shift the 
supply curve parallel to itself while a change in equilibrium price as a 
result of changes in d will pivot supply curve inwards or outwards 
depending on whether such a change is an increase or decrease in d. 

2. A certain consumer consumes two goods rice and beans. The utility 
derived from the two goods are given as follows: 

U = (3B 2 +BR+2R 3 ) 4 .3.8 
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Assuming that B = 3 and R = 2, determine the marginal Utilities of beans 
and rice consumed by the consumer. 

The rule to apply for finding the marginal Utilities of the consumer is rule 
3 (the generalized function rule). Using equations 3.6 and 3.7, 

5U/5B = 4(3B 2 +BR+2R 3 ) 3 (6B+R) 

¿U/oB = (24B+4R)( 3B 2 +BR+2R 3 ) 3 

SU/SR = 4(3B 2 +BR+2R 3 ) 3 (R+6R 2 ) 
oU/oR = (4R+24R 2 )( 3B 2 +BR+2R 3 ) 3 
With B = 3 and R = 2, 

SU/dB = (24(3)+4(2)][3(3) 2 +(3)(2)+2(2) 3 ] 3 
= [72+8][27+6+16] 3 
= (80)(49) 3 
= (80)(117649) 
oU/oB = 9,411,920 utils 

5U/5R = [4(2)+24(2) 2 ][3(3) 2 +(3)(2)+2(2) 3 ] 3 
= (8+96)(27+6+4) 3 
= (104)(37) 3 
= (104)(50653) 

¿U/oR = 5,267,912 utils 

3. Let the production function of a particular firm exhibit Cobb-Douglas 
production function of the form Q = AK“L^; (a+(3) 1. Find the 

marginal producís of labour and capital. The appropriate rule for 
solving this problem is the product rule. So we use equations 3.2 and 
3.3. 

MP k = — = (AZ/) — 

a k dk 

= AL p a K a_1 
= aAL p K “ 4 
_ aAL p K a 
K 
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MP k = 



MP l = ^=- = (AK a ) — Lf 
a L dL 

= (AK a )(3L p_l 

= AK a (3L p4 

= (3AKL P_1 

_ pAK a L p 

L 

MP k = J3- 
Z/ 

4. Let the utility function of a consumer exhibit constant elasticity of 
substitution (CES) utility function of the form: 

U = U[aXf p + (l-a)X 2 ' p ]- 1/p .3.9 

Where Xi and X 2 are commodity bundles consumed and a and [3 are 
structural parameters. Find the marginal utility of goods Xi and X 2 if a 
and (3 are 0.4 and 0.6 respectively and X¡ and X 2 are 10 and 8 units 
respectively. 

This problem requires invoking rule 3. Thus, we apply equations 3.6 
and 3.7 

-A = _a|3X 1 - ||W> 

ax í 

= "|3(l-a)Xo- (p+1) 

dx 2 

— = -l/(3[oX 1 - p + (l-a)X 2 ' p ] ( ' 1/p) ' 1 (-apXi- p ' 1 ) 
dx l 

= -1/ (3[aXf p + (l-a)X 2 ' p ]' (1+P)/p (-a pXi‘ (p+1) ) 

du/dx 1= aX!- (p+1) [oXE 15 + (l-a)X 2 ' p ]' (1+p)/p 
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du/dx 2 = -1/ p[aXf p + (l-a)X 2 " p ] ( - 1/P) ' 1 [|3(l-a)X2- (p+1) ] 

= -1/ (3[aXf p + (l-a)X 2 - p ]' (1+p)/p [|3(l-a)X2' (p+1) l 

du/dx 2 = [(l-a)X 2 ' (p+1) ][aX 1 ' p + (l-a)X 2 ' p ]- (1+p)/p 
With a = 0.4; (3 = 0.6 and X] = 10 while X 2 =8, then 
Marginal utility of good X] will be: 

du/dx 1= 0.4(10)' <0 ' 6+1) [0.4(10)'°' 6 + (l-04)10 OÍ T (1+0 ‘ 6)/0 ' 6 
= 0.4(10) 16 [0.4(10) °- 6 + (Í-O^IO' 0 - 6 ]- 1 - 670 - 6 

While marginal utility of good X 2 will be: 

du/dx 2 = [(l-0.4)(8)- ,06+1) ][0.4(10) 06 + (l-0.4)10- 0 - 6 ]- (1+0 ‘ 6)/0 ' 6 
= 0.6(8) 16 [0.4(10) °- 6 + (l-O^ílO 06 ] 1 - 670 - 6 

Total Differentials 

In our partial differentials, we looked at how a change in one parameter 
affects the endogenous variable, holding other parameters constant. In this 
section, we relax the assumption of holding other parameters constant. 
That is, we ask how the endogenous variable changes will give an 
infinitesimal change in all the parameters. 

Consider a saving function that is expressed as follows: 

S = S(Y,i).Ú.3.10 

Where S is total national savings Y equals aggregate income and i is the 
real interest rate. How will S change given an infinitesimal change in 
income and interest rate? This is solved for by applying the following 
formula: 

dS = (<9S/9Y)dY + (5S/5i)di 

dS = SydY + S¡di.3.11 

The expression dS, being the sum of the approximate change from both 
income and interest rate is called total differential of the saving function. 
The process of finding it is called total differentiation. We see that the 
sum of partial differential is actually the total differential of a function. In 
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a more general case, a function with independent variables can be 
exemplified by say, a utility function as in: 


U = U(X ls x 2 , X 3 ,., X n ) 


, TT a u , v a u , v du , v a u , v 

dU — -dX i + -dX 2 + -dX3 +.+ -dX n 

a.vj a.v 2 c)x 3 dx n 

=Xi- —— 1 dXi 
dx,. 


Rules of Total Differential 

Rule 1. dk = 0 

Rule 2. d(CU n ) = CnU n l dU 

Rule 3. d(U±V) = dU+dV 

Rule 4. d(UV) = Vdu + UdV 

„ . _ ,U VdU-UdV 
Rule 5. d — = 


V 




Examples 

2 

1. Find the differential dy of the function y = 5Xf + 3X 2 
We invoke rule 3. 


dy = 


dy 


dX í + 


dy 

\dX 2 j 


dX, 


.3.13 


dy = lOXidXi + 3dX 2 


(x H - X ) 

2. Find the total differential of the function y = —-— 

J r\ l 


2x: 


We make use of rule 5. 

a_y _ (Xj + 2X 2 ) 


We know that 


dX , 


2X; 


3.12 
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dy _ 1 

dX 2 ~ 2X\ 

ÍX,+2X 2 ) 1 

Henee dy = - — -^-dX ,+ -dX ? 

2X, 3 1 2X~ 2 


Total Derivatives 

In some cases, it is possible to discover that one of the independent 
variables is actually dependent on another independent variable in the 
model. For instance, it will be misleading to find a rate of change in 
consumption (C) with respect to tax, while holding income (Y) constant, 
because tax and income are, in some sense, related. A tool used to solve 
this problem is the use of total derivative. 

To make things easy, consider this function: y = f(X,W); X = h(W). 
We can combine the two functions together to get y = f(h(W),W) 

Indicating that the three variables are related to one another. W affeets y 
through direct and indirect channels. The indirect channel is through 
function h and the direct channel is through function f. thus: 

d^(dy_YdX_YdK , 1¿1 

dW dW . 



Application to Economics 

Let the production function of a particular firm be expressed as follows 
Q = Q(K,L,t) 
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That is, level of output depends on capital (K), labour (L) and time (t). 
But capital and labour can also depend on time as well (that is, they vary 
with time) such that 
K =K(t) L = L(t) 

The rate of change of output with respect to time can be expressed as 
follows: 


dQ 

dt 


'dQ\dK | 
y dK y dt 


'dQ\dL | 3 Q 
v 3 L ) dt dt 


3.15 


Example: let Q = A(t)KL = 4KL. K = 5+0.2t, L = 3+0.6t. Find the rate 
of change of 

output with respect to time 


c . , dQ dQ dK dL 3 Q 

First find-,-,- , — and -. 

dK dL dt dt dt 

dQ ,, dQ dK dL dQ 

— = 4L, — = 4K, -= 0.2,— = 0.6 and — = 0 

dK dL dt dt W dt 

Plug these in equation 3.15 

= 4L(0.2) + 4K(0.6) + 0 
dt 


dQ/dt = 0.8L+2.4K 

if the valúes of K and L are given, we can easily solve for the rate of 
change. 
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LECTURE FOUR 


Optimization Problems 


Introduction 

The common criterion of choice among alternatives in economics is for 
maximizing or minimizing problems. Conventionally, the primary 
objective of a firm is to maximize profit, while the basic objective of 
consumer is to maximize utility from consuming one or some 
combinations of commodity bundles. Also, firms are out to minimize cost 
of producing a certain level of output while consumers are willing to 
minimize expenditure on attaining given level of utility. The question 
that follows is: at what level of output are profits maximized or what is the 
minimum cost that can produce a specific level of output? In the case of 
consumers, we are faced with the questions of what is the level 
consumption combination for which máximum utility is realized or what is 
the minimum expenditure that will achieve a certain level of utility? 
When we are trying to find answer these questions, we are actually 
optimizing - seeking the best alternative. 

Broadly speaking, there are two categories of optimization - 
unconstrained optimization and constrained optimization. In this lecture, 
we shall be dealing with both constrained and unconstrained optimization 
and the methods of finding the máximum or minimum level of the choice 
variable. 

Objectives 

At the end of this lecture, you should be able to: 

1. identify conditions for which a function is minimized or 
maximized; 
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2. use the various methods for finding mínimum or máximum valué 
of a function; and 

3. solve some economic problems that has to do with optimization 

Pre-Test 

1. State the condition for which a function is maximized or 

3 2 

minimized. Given the function y = X + X' -3X. What is the level 
of X for which the function is maximized or minimized? 

2. What do you understand by constrained optimization? Maximize 
the foliowing production function Q = 3KL subject to 2 Pk + 5Pl = 
100 

CONTENT 

In optimizing any function, we need to identify the objective function and 
the choice variables. The objective function is the function that we want to 
maximize or minimize. The choice variables are the variable to be 
manipulated in order to get the expected result (either maximization or 
minimizadon). Consider the following profit function: 

n(Q) = R(Q) - C(Q).4.1 

Where Q is the output level, R and C are total revenue and total cost, 
respectively. Assuming we want to maximize equation 4.1 (the profit 
function). In this case, n is the dependent variable which represents the 
object to maximize, while R(Q) and C(Q) are independent variables or 
choice variables. Generally, given the function 

y = f(x).4.2 

We may attempt to develop a procedure for finding the level of x that will 
maxi mi ze or minimize the valué of y. 


JJnconstrained Optimization 
Máximum and Mínimum Valué 

We need to differentiate between absolute (or global) and relative (or 
local) máximum and mínimum valué. Consider the three graphs below. 
In panel A, the valué of x will result in the same valué of y, henee, the 
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height of each point A, B and C can be considered as máximum, mínimum 
or neither. This occurs when the function is constant. 

In panel B, we see that the function keeps on increasing. Point D is 
called absolute mínimum because there is no other point lower than that. 
However, since the function keeps on increasing, we cannot identify the 
máximum point. In panel C, we cannot consider point E as the global 
mínimum. Also, we cannot consider point F as the global máximum 
because the function may either decrease or increase further. Thus, points 
like E and G are called relative (local) mínimum, while point like F is 
called relative (local) máximum. In most economic problems, we 
usually deal with relative extremum of this type. 



Relative extremum is a point at which a function is at a relative máximum 
or mínimum. It must be added that a point in the domain of a function 
where the derivative equals zero or is undefined is called a critical point or 
valué. 

We use the concept of derivative to determine whether a function is at 
its relative mínimum or máximum. For convenience, we shall refer to the 
derivative of function as first (order) derivative. For an extremum to exist, 
the necessary condition is that the first derivative must be zero. The valué 
of the function when its first derivative is set to zero is called the critical 
point. A relative extremum must also be a stationary valué, but a 
stationary valué may be associated with either a relative extremum or 
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inflection point. Stationary valué is the valué that a function assumes at its 
critical valué. An inflection point occurs at the point where no particular 
direction of a slope is determined. 


Unconstrained Optimization with univariate function 

To find the relative máximum or minimum of a univariate function, (f(x), 
take these steps: 

a. take the first order derivative and set it equal to zero (first order 
condition) 

b. obtain the critical point by evaluating the valué of f (x) = 0 

c. you may find the stationary valúes of the function. This is 
established by plugging the value(s) of the critical point in f(x) 

d. take the second derivative, f"(x) (second order condition) and 

e. evalúate it at the critical point(s) and check the sign(s). If at a 
(particular) critical point x, f"(x) < 0, the function is at a relative 
máximum (concave). If f"(x) > 0, the function is at a relative 
minimum (convex). 

The first order derivative gives the necessary condition while the second 
order derivative gives the sufficient condition. 

Example 1 ^ V 

Optimize: f(x) = 2X 3 -30X 2 +126X+59 

Take the first order condition to establish the necessary condition 
f(x) = 6X 2 -60X+126 = 0 

Evalúate this function to get the critical point. The function has turned to 
quadratic equation and thus to evalúate, we shall use quadratic formula. 
This can be done by factorization or by adopting the general formula for 
quadratic equation. In this example, we use the former: 

6(X-3)(x-7) = 0 
X = 3, X = 7 (critical points) 

Take the second order derivative, evalúate it at the critical points and 
determine the status of the function. 
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f"(x) = 12X-60 

f"(3) = 12(3)-60 = 36-60 = -24 < 0 (relative máximum) 

f"(7) = 12(7)-60 = 84-60 = 14 > 0 (relative minimum) 
to obtain the stationary point, plug X = 3 and X = 7 in the f(x) function. 

We conclude, therefore, that the function (fx) is at máximum when X is 3 
and it is at minimum when X is 7. 

Example 2 

Find the relative extremum of the function: 
y = f(x) = X 3 -12X 2 +36X+8 
f'(x) = 3X 2 -24X+36 = 0 
= X 2 -8X+12 = 0 

factorize 

X 2 -2X-6X+12 = 0 
X(X-2)-6(X-2) = 0 
(X-6)(X-2) = 0 

Henee X = 6, X = 2: (critical points) 
f"(x) = 6X-24 

f"(6) = 6(6)-24 = 36-24 - 12 > 0 (relative minimum) 

fff( 2 )V 6(2)-24 = 12-24 = -12 < 0 (relative máximum) 
we thus conclude that the function f(x) is at relative minimum when X=6 
and at máximum when X=2. 

Application to Economics 

1 Find the optimal valué for average cost when: 

AC = f(Q) = Q 2 -5Q+8 

Where AC is Average cost and Q is level of output. 
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Take the first order derivative (FOD) = the necessary condition 
f (Q) = 2Q-5 = 0 
Q = 5/2 = 2.5 (critical point) 


The second order derivative (20D) = the sufficient condition 


f"(X) = 2 > 0 (relative minimum) 

henee the minimum average cost (MAC) is attained when output is 2.5 
units 

2 Maximize the profits n for a firm, given total revenue and total cost 
functions: 

R = 4000Q-33Q 2 

C = 2Q 3 -3Q 2 +400Q+5000 

Where R and C are total revenue and total cost, respectively, and Q is 
level of output. 

Determine the profit function first using the formula n = R(Q)-C(Q) 
n = 4000Q-33Q 2 -[2Q 3 -3Q 2 +400Q+50001 
= -2Q 3 -30Q 2 +3600Q-500 
Now take FOD of the profit function 
n = -6Q 2 -60Q+3600 = 0 
= -Q 2 -10Q+600 = 0 

Applying general formula for quadratic equation, we find out that the 
critical points are 

Q = -30, Q = 20. 

We need to determine which of the Qs actually maximizes profit. This is 
done by taking the 20D 


n = -12Q-60 

n (-30) = -12(-30)-60 = 360-60 = 300 > 0 (relative minimum) 
n (20) = -12(20)-60 = -240-60 = -300 < 0 (relative máximum) 
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henee, the máximum profit is attained when output level reaches 20 units. 
What is the profit at this level of output? To find this, substitute 20 for Q 
in the profit function 

n(20) = -2(20) 3 -30(20) 2 +3600(20)-500 

= 39000 

Meaning that the optimal profit is 39000 and the level of output that 
generates this profit is 20 units. 

3. a firm has the foliowing total cost and demand functions: 

C = ^ Q 3 -7 Q 2 +111Q+50 
Q d = 100-P 

Where C and Qd are total cost and quantity demanded and P is unit price 
of output Q. 

Formúlate the total profit function. Determine the profit maximizing 
level of output Q . What is the máximum profit at this level of output? 

Formúlate the profit function as follows 
n = R(Q) - C(Q) 

R(Q) = PQ 

Thus n = PQ-C(Q) 

From the problem given, P = 100-Q (why?) 

Thus R(Q) = (100-Q)Q = 100Q-Q 2 

n = 100Q-Q 2 - (^Q 3 -7Q 2 +111Q+50) 

= 100Q-Q 2 - ^ Q 3 +7Q 2 -111Q-50 
= - ~ Q 3 -Q 2 +7 Q 2 +100Q-111 Q-50 

n = - ^ Q 3 +6Q 2 -llQ-50 = profit function 
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To find the profit maximizing level of output, take the FOD of the profit 
function and find the critical point 

n = -Q 2 +12Q-11 = 0 

Simplifying using general formula for quadratic equation, we have 

Q = 1, Q = 11 (critical points) 

Take the 20D to establish which of the two valúes of Q satisfies sufficient 
condition for máximum 

n = -2Q+12 

n (1) = -2(1)+12 = -3+12 = 9 > 0 (relative minimum) 

n (11) = -2(11)+12 = -22+12 = -10 < 0 (relative máximum) 

Thus, the profit maximizing level of output is 11 units 

to find the máximum profit at this level of output, insert 11 into the profit 
function 

n = -1/3(11) 3 +6(11) 2 -11(11)-50 
= -1/3(1331)+726-121 -50 
= -443(l/6)+726-121-50 

n = 140 or 111- 

. < 3 

Unconstrained Optimization with Multivariate functions 

We may encounter a situation whereby the objective variable is a function 
of more than one decisión variables. In this case, we need to modify our 
rules in order to capture this scenario. 
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For a multivariate function such as z = (x,y) to be a relative mínimum 
or máximum, three conditions must be satisfied and these are: 

a. the first-order pardal derivative (FOPD) must equal to zero 
simultaneously 

b. the second-order pardal derivative (20PD), when evaluated at the 
critical point must both be negative for a relative máximum and 
positive for a relative mínimum 

c. the product of the 20PD evaluated at the critical point(s) must 
exceed the square of the second-order cross-partial derivative 
(20CPD) also evaluated at the critical point. 

The importance of the third condition is to preelude the inflection or 
saddle point. In summary, for a 2-variable function, say z = (x, y), 


f x , f y = 0 (necessary condition) 
f xx , fyy < 0 condition for máximum 
f xx , fyy. > 0 condition for mínimum 

fxx fyy > (fxy) sufficient condition for relative máximum or mínimum 


Examples 

o T 

Optimize the following unconstrained function: z = 2y -x +147x-54y+12 
Take the FOPD and set them at zero 

z x = -3 x 2 +147 = 0 z y = 6y 2 -54 = 0 

henee 3x 2 = 147 6y 2 =54 



x = ± 7 y = ± 3 

We now have four (4) critical points viz: 
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Next, take the 20PD, evalúate them at each critical point and check the 
relative extremum: 

z xx = -6x z yy =12y 

Z x y 0 Zyx = 0 

Z xx (7,3) =-6(7) = -42 < 0 z yy (7,3) = 12(3) = 36 > 0 

Zxx (-7,3) = -6(-7) = 42 > 0 z yy (-7, 3) = 12(3) = 36 > 0 

Zxx (7,-3) = -6(7) = -42 < 0 z yy (7,-3) = 12(-3) = -36 < 0 

Zxx (-7,-3) = -6(-7) = 42 > 0 z yy (-7,-3) = 12(-3) = -36 < 0 

Now take the third and sufficient condition and verify which of the 20PDs 
meet the condition: 

-42.36 < (0) violates the rule thus we strike this out. 

42.36 > (0)“ satisfies third order condition. The critical points are 
a candidate 

-42.-36 > (0)“ satisfies third order condition. The critical points are 
a candidate. 

42.-36 < (0) violates the rule thus we strike this out. 

Two sets of the four critical points satisfy the sufficient condition. The 
first one is when x = -7 and y = 3. The second is when x = 7 and y = -7. 

This implies that the function is at mínimum when x = -7 and y = 3. 
The function is at máximum when x = 7 and y = -3. 

Constrained Optimization 

Up to this point, we assume that the objective function being optimized 
has nd constraints. In real life, particularly in economic analysis, there are 
several constraints that face a given objective. For instance, the optimal 
level of utility derived from a particular set of commodity depends on the 
income level of the consumer and the prices of the goods among others. 
Thus, income and prices of goods serve as constraints to consumer’s 
satisfaction. We therefore need to recognize these limiting factors and 
incorpórate them in our optimization problem. 
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Consider a consumer with a simple utility function of the form: 

U = XiX 2 + 2Xi.4.1 

It is straightforward to verify that without any constraint, the consumer 
will tend to continué consuming the goods until he is tired. Assuming that 
the utility function is subject to the purchasing power of the consumer (the 
budget constraint). Assume further that the consumer intends to spend 60 
naira on the two goods and the current prices of the goods are 4 naira and 
2 naira, respectively. The budget constraint is written as follows: 


4Xi+2X 2 = 60.. 4.3 

What will be the commodity bundles that maximizes consumer’s utility? 
We can solve for this using three methods: 

• substitution method; 

• LaGragian-multiplier method; and 

• Bordered Hessian method. 


The Use of Substitution Method 

Solve for X¡ in terms of X 2 in equation 4.2 and substitute this in equation 
4.1 


x _ 60-2X 2 


4.4 


Henee, U = X 2 


60- 


■2X 1 60-2X^ 

2 + 2 - 


4 


4 


Simplifyinf', we have 56X 2 - 2X 2 +120.4.5 

Now take the FOD and set it equal to zero to obtain the critical point 

U x2 = 56-4X 2 = 0.4.6 

The critical point gives 


X 2 =56/4 = 14 

plug this valué in 4.3 to obtain the critical valué of Xi 
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x = (60-2)4 _ 32 _ g 
1 4 4 

The stationary valué can be found by inserting the critical points into the 
objective function (equation 4.1) 

U* = (8) (14)+2(8) = 112+16 = 128 
To determine the relative extremum status, take the 20D: 

Ux 2 X 2 = -4 < 0 (relative máximum) 

Thus, given the budget constraint, the máximum level of utility is attained 
when the consumer buys and consumes 8 units of Xi and 14 units of X 2 . 
The máximum utility attainable from the consumption of these good is 

128 utils. 

Example 2 

Let the production function of a firm be Q = 4K“+3KL+6L“. let the per 
unit costs of L and K be 1 naira each and the amount available for the 
purchase of these factor inputs is 56 naira. Determine the optimal level of 
output. 

Set up the objective and cost constraint problems 


Max 4K 2 +3KL+6L 2 ..r. .4.7 

St K+L = 56.4.8 

By substitution method, using equation 4.8, 

K = 56-L.4.9 

Substituting equation 4.9 in 4.7, we have 


Q = 4(56-L) 2 +3(56-L) + 6L 2 

112664-448L+4L 2 +156L+6L 2 


12664-282L+10L 2 .4.10 

Now take the FOD of equation 4.9 

Ql = -282+20L = 0.4.11 
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4.12 


L = 282/20 14.1 = 14 . 

Substitute equation 4.12 for L in equation 4.9 

K = 56-14.1 = 41.9 = 42 .4.13 

To find the stationarity point, plug equations 4.12 and 4.13 in 4.7 
Q = 4(42) 2 +(3) (42)(14)+14 2 
= 7056+1764+196 
Q = 9016 

Thus, given the cost constraints faced by the firm, the level of input 
combinations that maximizes output are 42 units of capital and 14 units of 
labour and the optimal level that can be produced is 9,016 units. 

The use of LaGragian Method 

We can handle the constrained optimization with the use of LaGragian 
multiplier. 

Example: Use LaGrangian multipliers to optimize the following functions 
subject to the given constraint: 

Z = 4X 2 -2XY+6Y 2 subject to the constraint X+Y = 72 

Steps to take 

• Set up the constraint equal to zero, multiply it by X, and add it to 
the objective function 

• Take the FODs of the function and set it to 0 

• Obtain the critical points 

• Put the valúes of the critical points into the primitive function to 
obtain the stationary point 

• Test for the necessary condition using the same process discussed 
earlier 

Taking the first step, we have: 


Z = 4X 2 -2XY+6Y 2 + M72-X-Y) .4.14 

Second step: 

Z x = 8X-2Y-7, = 0.4.15 

Z Y = -2X+12Y-7, = 0.4.16 
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= 72-X-Y=0.4.17 

Third step: 

From equations 4.15and 4.16, 

8X-2Y = -2X+12Y 
Henee 10X = 14Y 

Y = 10/14 = 0.7X.4.18 

Substitutive equation 4.18 into 4.17 
72 = X+0.7X 
72= 1.7X 

X* = 42.5 .4.19 


Substitute 4.19 into equation 4.18 to obtain the valué of Y 

Y* = 1.7(42.5) = 29.75 

Substitute the valúes of X and Y into equation 4.16 to obtain the 
stationary point 

Z = 9936. 

Utility Maximization with Constraint 

Consider the utility of a consumer described as follows: 

U = U(X,Y); (U X ,U Y > 0) 

Subject to P xX+P y Y = I 

Maximize the utility subject to income constraint 


> Set up the LaGragian function 

L = U(X, Y) + X (I- PxX-PyY).4.20 

> The FOD of 4.18 yield the following 

L x = U X -AP X = 0.4.21 

L y = U Y -AP Y = 0.4.22 

U = I- P x X-P y Y = 0.4.23 

Equations 4.19 and 4.20 are equivalent to 
U P 

.4.24 

U Y P Y 
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4.25 


Rearranging we have 

u JL= u^ =l 

P X Py 

Equation 4.23 says that in order to maximize utility, consumers must 
allocate their budgets so as to equalize the ratio of marginal utility to own 
price for every commodity. The term X is the marginal utility of money. 

Profit Maximization and Cosí Constraint 
Let the profit function of a firm be of the form 
II = fI(K,L) subject to wL+rK = Y 

Where fl is the profit, K and L are quantities of capital and labour while w 
and r are wages rental rates. The máximum profit is attained as follows: 

Z = n (K, L) + X(Y-wL-rK).<.4.26 

The FOD of the function gives: 

Zk = nK- 7,r = 0.,.4.27 

Z L = nL-^w = 0...4.28 

Z x = Y-wL-rK = 0.r^. % .4.29 

From equations 4.27 and 4.28, we generate equation 4.30 below 

n,„ /- 

— = - = X . 4.30 

n L w 

The interpretation of equation 4.28 is similar to utility maximization 
problem discussed earlier. 

Using Bordered Hessian to handle Optimization Problems 
Elementary understanding of Bordered Hessian 

In optimization, most especially in the case of multivariate equations, 
there are 3 conditions to be met. Given =y = (x¡, X 2 ) 

First condition: f x i= 0; f X 2 = 0 

Second condition: f x i x i, f X 2 x 2 > 0 relative minimum 
fxixi, fx 2 x 2 < 0 relative máximum 
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Third condition: f x i x i.f x2x2 > (f x ix 2 ) 2 


One simple way for establishing the second condition, particularly when 
we are dealing with more than two variables is the use of bordered 
Hessian. A Hessian is a determinant composed of all the second order 
derivatives. While the second order direct partíais are on the principal 
diagonal, the second-order cross partials are on the off the principal 
diagonal. 

To make things easy, let us, for a moment, consider a bivariate 
function of the form y = f(xi,X 2 ). The second order derivative and the 
second order cross partial derivative can be arranged in matrix format as 
follows: 


f 

V 


fx ni 

fx2x\ 


fx\x2 

/, 


x2x2 J 


Henee IHI = 


fX lxl 

fx 2x1 


fx\x2 


x2x2 J 


We need to solve for the determinants of the 
matrix above, there are two principal minors. 
deleting the last row and the last column. This 
the principal diagonal, that is: 


principal minors. In the 
The first one is got by 
gives the first element on 


IHd = f xbd 


The second principal minor is given by 

( f f > 

J x\x\ J x\x2 
\fx2x\ fx'. 


ih 2 i 


' x2x2J 


- fxlxl.fx2x2- fx 1 X 2- f X 2 X l 
= fxlxl.fx2x2-( f X l X 2) 2 
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Rules for máximum and mínimum for bivariate function 
If the valué of IHjl and IH 2 I are positive then, the Hessian IHI is called 
positive definite. A positive definite Hessian satisfies the second-order 
condition for relative minimum. 


If the valué of IHil is negative while the valué of IH 2 I is positive, then 
the Hessian IHI is called negative definite. A negative definite Hessian 
satisfies the second-order condition for relative máximum. 

Let us now extend our analysis to three variables. Let y = f(xi,X 2 ,X 3 ) 
The second-order for the variables is: 


fxlxl = 

yn 

fx2x2 = 

yn 

fx3x3 = 

= y33 

fxlx2 = 

: y2i 

fxlx3 = 

Y31 

fx2x3 = 

= y32 

fx2xl = 

: yi2 

fx3xl= 

yi3 

fx3x2 = 

= yo 


hn 

y 12 

A 3 ' 



IHI = 


y 22 

y 23 




lAl 

y 32 

y 3 3 j 




Next, we need to find the determinants of the principal minors. Since this 
is a 3x3 matrix, we shall have 3 principal minors. The first is got by 
deleting the last two columns and the last two rows. 

IHil = yn 

To get the determinant of the second principal minor IH 2 I, we delete the 
last column and the last row of equation 4.31 



Ai y 12 

_y 2 \ y22_ 

yny22-(y2i) 2 


The IHI is the third principal minor 
IH 3 I = IHI. 
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Rules for mínimum and máximum for multivariate function 

If the valué of IHi I, IFW and IH3I are greater than zero, then the Hessian is 

positive definite and henee satisfies condition for minimum. 

If IHi I < 0 , IH2I > 0 and IH3I < 0 , then the Hessian are negative definite 
and henee satisfies condition for minimum. 

Bordered Hessian and Unconstrained optimization 

Let the utility function of a consumer be represented by this behavioural 
function: 

U = -5Xi 2 +10Xi+XiX 3 - 2Xj 44 X 2 + 2 X 2 X 3 - 4X 2 

Determine the level of consumption mix that maximizes the consumer’s 
utility. 

The first-order conditions are: 

Ui = -10Xi+10+X 3 = 0 
U 2 = -4X 2 + 2 X 3+4 = 0 
U 3 = Xi+2X 2 -8X 3 = 0 

-10Xi+ X 3 = -10 
-4X 2 +2X 3 = -4 
Xi+2X 2 -8X 3 = 0 

This can be written in matrix form as in: 


"-10 

0 

n 

(xA 


'- 10 ' 

0 

-4 

2 

x 2 

= 

-4 

V 1 

2 

- 8 , 

UJ 


l 0 J 


With the use of matrix inversión or Cramer’s rule, we find out that 
X\ = 1.04; X* 2 = 1.22 and X* 3 = 0.43 
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To be sure that these commodity bundles satisfy the maximization 
condition, we need to test for the second order condition. This can be 
done with the use of Bordered Hessian: 

Un = -10, U 12 = 0 Ui 3 =l 
U 2 i = 0 U 22 = -4 U 32 = 2 

U 31 = 1 U 32 = 2 U 33 = -8 

'-10 0 1 
IHI =0-42 
v 1 2-8 

IHil = -10 < 0 




= -4(-10)-0 = 40 > 0 
|H 3 I = IHI = -10(28)44 = -276 < 0 

Thus, the utility is maximized at the consumption level of X*i = 1.04; X *2 
= 1.22 and X * 3 = 0.43. 

Bordered Hessian and Constrained Optimization 

Let the production of a firm exhibit Cobb-Douglas of the form: 

< Q = K° V 5 

Let the cost per unit of capital (K) that is Pk be 3 while the per unit cost of 
labour (L) that is Pl be 4 and the total expenditure is 108. Determine the 
level of input combinations that gives máximum output. 

Max K°V- 5 

St 3 P k + 4P l =108 
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Set up the LaGragian function: 


Z = k°- 4 L°' 5 +>u( 108-3P k - 4P l ) 

Z k = 0.4K"°' 6 L°' 5 -3A, = 0 

Z L = 0.5K°' 4 L"°' 5 -4A, = 0 

Zx= 108-3 Pk-4Pl=0 

Solving for critical points, we have: 

K = 16, L = 15 

To test for the second-order condition, we use Bordered Hessian. From 
the above information, 

Zkk = -0.24K~ 16 L 05 
Z ll = -0.25K°V- 5 
Z LK = Zkl = 0.2K"°' 6 L " 0 ' 5 
Zpk = 3 
Z PL = 4. 

Henee 

'-0.24 R- L6 L 05 O.IK^U 0 - 5 3 

IHI = 0.2K-°- 6 nÉ- -0.25K OA L~° 5 4 

3 4 0 

V 

IHjl = -0.24K" 1 ' 6 L 0 ' 5 < 0 

IH 2 I = (-0.24K" L6 L 0 ' 5 ) (-0.25K 0 ' 4 L 4 ' 5 )-(0.2K'°' 6 L' 0 ' 5 ) (0.2K'°- 6 L' a5 ) 

= 0.06K'°' 2 L-0.04K" 1 ' 2 
= 0.06L/K°' 2 -0.04/K L2 > 0 
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Solving for IH 3 I, we find out that the valué is actually < 0, thus, we have 
IHi I < 0, IH 2 I > 0 and IH 3 I < 0. We, therefore, can conclude that when the 
capital is hired up to 16 units and the labour is hired up to 15 units, 
máximum output will be realized. 


59 



LECTURE FIVE 


Function Approximation 


Introduction 

Up to this point, we have not considered a situation whereby the second 
derivative turns out to have zero valué at the stationary point. If this is 
the case, none of the methods discussed earlier could handle the situation. 
This chapter will discuss how such problem can be handled. To be able to 
handle this task appropriately, we need to intimate ourselves with 
expansión of a function around point x = 0 or around any point x = xo 
where xo is a specific valué not equal to zero. The former is called 
Macluarin series, while the latter is called Taylor series. 

To expand a function y = f(x) around zero (x=0) means to transform 
that function into a polynomial form, in which the coefficients of the 
variable terms are expressed in terms of the derivative valúes f(x), f (x), 
f"(x), f'"(x) etc, ah evaluated at the point of expansión, 0. For the Taylor 
series, the evaluation is done at the point of expansión xo, i.e. f'(xo). f"(xo), 
f'■"(xo) etc. The result of expansión is a power series because being a 

polynomial; it consists of a sum of power functions. 

Function Approximation of Maclaurin Series 

Consider the polynomial function of the nth degree 

f(x) = ao+ aix+a 2 X”+a 3 X +a 4 X +.+ a n x .5.1 

the power series after expansión will involve the derivation of various 
orders of the function f. Thus, by successive differentiation of equation 
5.1, we have: 
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f (x) = ai+2a2X+3a3X 2 +4a4X 3 + ... 

.+na n x nl . 

. 5.2 

f"(x) = 2a 2 +3(2)a 3 x+4(3)a4X 2 + . 

.+n(n-1 )a n x n 2 .... 

..5.3 

f "(x) = 3(2)a 3 +4(3)(2)a 4 x+. 

.+n(n-l)(n-2)a n x n ‘ 3 ... 

. 5.4 

fVx) = 4(3)(2)a 4 x+. 

. +n(n-l)(n-2)(n-3)a n x n ~ 4 . 

.5.5 

f 4 (x) = (n-l)(n-2)(n-3). 

.3(2)(l)a n . 

... 5.6 


We discover that successive differentiation reduces the number of 
terms by one until, in the nth derivative; we have a single product term 
which is constant. These derivatives can be evaluated around any valué, 
but the case we are considering is its evaluation around x V 0. For x = 0, 
all the constant terms in equations 5.2 through 5.6 will drop out and what 


will be left is 

f (0) = ai. x ....!.5.7 

f"(0) = 2a 2 ...5.8 

f "(0) = 3(2)a 3 ..y;...5.9 

f 4 (0) = (4)(3)(2)a 4 .v;...5.10 


f n (0) = (n-l)(n-2)(n-3).(3)(2)(l)a n .5.11 


Let us adopt a shorthand to reduce this long term. Assuming that 

n! (Pronounce it as n factorial) = n(n-l)(n-2)(n-3).(3)(2)(1) such that 

2! will be 2x1 and 4! will be 4x3x2xl = 24. A special case is where 0! = 1. 
Note that we can write equations 5.7-5.10 as 

f (0) = a, 

/"( 0 ) _ 


2 ! 

/'"(0) 
3(2)! 

/ 4 ( 0 ) 
4(3)(2)(1)! 


&3 


a 4 
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_ f\ 0 )_ 

n(n-l)(n-2)(n-3).(3)(2)(1)! 


thus aj= .5.11 

1! 

a 2 =m.5.12 

2 ! ^ 

a 3 =^. 

3! 

a 4 =/^ . jé*.. .5.14 

4! .vjs 

a n =m . M .5.15 

n ! 

if we substitute equations 5.11 through 5.15 into 5.1, we ha ve: 

/'(O) /"(O) 2 /'"(O) 3 f\ 0) 4 /"(O) „ 

f(x)=f(0)+ x + x + -——¿T + x + .+ X ..5.16 
1! 2! 3! CV 4! n! 


This new polynomial function f(x), which represents the expansión of 
the function f(x) around zero (x=0), is called Maclaurin series of the 
polynomial function. The important thing to note is that the point of 
expansión (0) is the valué of x that will be used to evalúate f(x) and all its 
derivatives. 


Examples 

2 

1. Find the Maclaurin series for the function: y = f(x) = 2+4x+3x 
First take the derivatives: 
f (x) = 4+6x henee f (0) = 4 
f"(x) = 6 henee f"(0) = 6 
Thus, the Maclaurin series is: 


f(x) = f(0)+f(0)x+^p.x 2 


= 2+4x+ — x 2 
2 ! 


= 2+4x+3x 2 
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The last two lines are used to verify that the Maclaurin series correctly 
represents the function. 


2. Given the function: y = h(x) = - = (1+x)" 1 , find the Maclaurin 

1 + x 


series 


h'(x) = -( 1 +x) 2 henee h'( 0 ) = -1 
h"(x) = 2 (l+x )' 3 henee h"( 0 ) = 2 
h"'(x) = -óü+x )" 4 henee h"'( 0 ) = -6 
h 4 (x) = 24(l+x)' 5 henee h 4 (0) = 24 


,, . , h\ 0 ) h "( 0 ) 2 

h(x) = h( 0 ) + -^^x +—— x 


1 ! 
2 _ 

1 ! 2 ! 


, -1 2 2 - 6 3 
— 1+ - H-X + -X + 


2 ! 

-6 

3! 


, h"\ 0 ) 3 , h\0) 4 

\ -x -l-x 

3! 4! 

24 4 
— x 
4! 


3. Pro ve that the Maclaurin series of the function g(x) = e x is equal to 
2.71828. 

Generate the derivative of various orders for the function. Note that 
the derivative of e x is e x , henee that the derivative of g(x) is simply e x . 
The second, third up to n derivative will be e x , henee, when we evalúate 
all the derivatives at the expansión point x= 0 , we have: 

g'( 0 ) = g”( 0 ) = g"'( 0 ) =.= g n ( 0 ) = e° = 1 . 

By setting xo = 0, the Maclaurin series of e x is 

e x = h(x) = h( 0 )+h'( 0 )x+h"( 0 )x 2 /2 !+h"'(0)x 3 /3!+.+h n ( 0 )/n! 

the Maclaurin series may then be expressed as a convergent infinite series 
as foliows: 

Y .. I 9 I 9 1 1 le 

& ~ l+XH- X~ H- X H- X H- X + . 

2! 3! 4! 5! 


for x= l, 

, . 1 1 1 1 

e = 1 + 1 + —+ — + — + —+. 

2! 3! 4! 5! 

= 2+0.5+0.1666+0.041666+0.008333+ 

= 2.7182 
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The important thing to note about this expansión is the appropriate rule of 
derivative we need to invoke and the number of times we need to 
differentiate. Also, note that rules of differentiation to adopt depend on 
the nature of the function. In other words, it is possible to invoke 2 or 3 
rules in an expansión. 

Function Approximation of Taylor Series 

As said earlier, it is possible to expand an equation at any point xo, not 
necessarily zero. Without loss of generality, let us use examplel as our 
guide in explaining Taylor series. For pedagogy, let us import the 
equation down here: 

y = f(x) = 2+4x+3x 2 
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For expansión around a specific point xo, let us assume that a given valué 
of x is actually the deviation of xo. In this case, x is equal to the deviation 
valué, say Ax (or 5) and xq, i.e. 


x = xo+Ax 


5.16 


If we substitute the RHS of 5.16 for x in our example, we obtain another 
function as follows: 

y = f(x) = 2+4(xo+Ax)+3(xo+Ax) 2 .5.17 

f (x) = = 4+6(x 0 +Ax) 

f"(x) = 6 A' 

Now, we know that Xo is a fixed number (a chosen number), henee, Ax is 
the only variable in the function. It follows that f(x) is a function of Ax. If 
we let h(x) = Ax then 

h(Ax) = 2+4(x 0 +Ax)+3(x 0 +Ax) 2 = f(x) 
its derivatives are: 

h'(Ax) = 4+6(x 0 +Ax) = f'(x) 
h"(Ax) = 6 = f"( x) 

Our understanding of expansión around zero indicates that the above 
expansión yields 


h(Ax) = 




- + - 


.5.18 


0 ! 1 ! 2 ! 
if Ax = 0, then, according to equation 5.16, x = xo. Therefore, on the basis 
of identity, 


f(0) = f(x 0 ), f (0) = f (x 0 ), f"(0) = f'(x 0 ) 

If we substitute this into equation 5.18, we have: 


h(Ax) = f(x) = 


f(x 0 ) + f^) f^) (x _ Xo)2 


0 ! 


1 ! 


2 ! 


For the function under consideration, 
f(x 0 ) = 2+4x 0 +3x 0 2 
f (x) = 4+6x 0 
f"(x) = 6 
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Therefore, applying equation 5.18, the Taylor polynomial expansión 
becomes 


f(x) = 2+4xo+3x 0 2 +(4+6xo)(x-xo)+6(x-xo) 2 /2 

2+4x+3x 2 


The generalization to the n-th degree polynomial of equation 5.18 gives 
the Taylor’s formula given in equation 5.19 


f(x) = 


f ( x o).f(x 0 ),_ , f'(x 0 ) 


0! 


- + - 


1! 


" (x _ Xq)+ - 


2 ! 


(x-x 0 )+. 


.+ 


f n (x 0 ) 


n ! 


(x~x 0 ) 5.19 


The difference between Maclaurin series and Taylor series is the 
replacement of zero by xo at the point of expansión and in the replacement 
of x by (x-xo). The interpretation of equation 5.19 is that given n-th 
degree polynomial f(x), if we let x = 5, say, in the terms on the right of 
equation 5.19, if we select an arbitrary number xo, then evalúate and add 
these terms, we will end up with exactly f(5) - the valué of f(x) at x = 5. 
Let us assume that xo = 4. we can write equation in example 1 as: 

fY4Y f "(41 

y = f( x ) = f(4)+f(4)(x-4)+ ——^-{x-A) 2 + .+ (x - 4)" 

2! ni 


Taylor’s Theorem 

This theorem says that given an arbitrary function p(x), if we know the 
valué of the function at x = xo (p(xo)), and the valúes of its derivatives at 
xo, p'(xo), p"(xoV p'"(x 0 ) etc., then this function can be expanded around 
the point xq as follows: 


M(x)-^+^(x-x 0 )+^(x-x 0 ) 2 +^(x-x 0 y+. 


0! 


1! 


2 ! 


3! 


ni 

.5.20 


- Pn+Rn 


Where P n represents the terms in the bracket and R n denotes a remainder. 
Because equation 5.20 bears R n , it is called Taylor’s formulation with 
remainder. The form of P n and the size of R n depend on the chosen valué 
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of n. The larger the n, the more terms there will be in P n , while R n will 
assume different valué for each different n. 

The appearance of R n is used to capture the polynomial expansión of 
arbitrary function p which cannot always be transformed exactly into, but 
can only be approximated by the polynomial form of Taylor’s formula. 
Thus, R n supplements for the discrepancy between p(x) and P n . 

If we choose n=l, we will have 

p(x) = p(x 0 )+ p'(x 0 )(x-x 0 )+Ri = P 1 +R 1 .5.21 

a linear approximation 
with n = 2 

p(x) = p(x 0 )+ p'(x 0 )(x-x 0 )+ p(x-x 0 ) 2 + R 2 = P 2 +R 2 .5.22 

a quadratic approximation 

Example 

Expand the function p(x) = 1 around the point xo = 4 

1+ X 

P (1) = 1/2 

p'(x) = -(1+x)' 2 henee pi'(l) = -(l)' 2 = -1/4 
p"(x) = 2(l+x)~ 3 henee p"(l) = 2(2) -3 = 14 
p'"(x) = -6(l+x)~ 4 henee p'"(l) = -6(2)' 4 = -3/8 
p 4 (X) = 24(l+x)’ 5 henee p 4 (l) = 24(2)' 5 = 3/4 

Setting xo = 1 and utilizing the obtained derivative, we get the following 
Taylor series with remainder 

p(x)=---(x-l) + -(x-l) 2 + —(x-l) 3 + —(x-l) 4 +R 4 
2 4 8 16 32 

Using Pascal triangle for polynomial expansión, we have 


, , 31 13 1 2 3 3 1 4 _ 

u(x) =-x +—x- X -t-X +R 4 

32 16 2 16 32 


if we know the valué of x, then we can get the approximate valué of the 
p(x) expanded at xq = 1 and n = 4 
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Exercises 


1. Expand the function p(x) = -- around the point Xo = 3 

l + x 

2 

2. Consider the function v(y) = (l-y')(y+4). Expand the point around 
x 0 = 4 


, 
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LECTURE SIX 


Integral Calculus (Indefinite Integral) 


Introduction 

This is the process of tracing back a given derivative to its primitive 
function. In other words, it is a way by which we move from the 
differentiated function to its original function that is, reversing the process 
of differentiation by finding the original or primitive function, from the 
derivative. This is also called integration or anti-differentiation. Thus, the 
original function f(x) is called the integral or anti-derivative of F 1 (x). 

For example, if the differentiation of a given primitive function F(x) 
gives the derivative f(x), then, we can “intégrate” f(x) to get the primitive 
function f(x) back, as long as we have the appropriate information 
available to definitize the arbitrary constant that will arise in the course of 
the integration. Invariably, integration involves the tracing of the 
parentage of the function f(x), whereas differentiation shows the progeny 
of the function EJx). At this point we must note this difference that, while 
the differentiation of the primitive function F(x), produces a lone 
offspring, namely, a unique derivative f(x), the derived function f(x) can 
be traced to an infinite number of possible parents through integration, 
because F(x) is an integral of f(x), then, also must be F(x) plus any 
constant. 

The standard notation to denote the integration of f(x) with respect to 

x is 

| f(x)dx .6.1 

The left side symbol that looks like an elongated S is called integral 
sign, while the f(x) part is known as the integrand (that is, the function to 
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be integrated), and the dx shows that the operation has to be performed 
with respect to the variable x, this integral sign that is in the front of the 
function f(x) dx gives the instruction to reverse the differentiation process 
that gave rise to the differential. We can write that 

J f(x)dx = F(x) + C . 6.2 


The C is the arbitrary constant of the integration, which indicates the 


múltiple parentage of the integrand. 


This integral 


jf(x)dx 


is the 


indefinite integral of f(x) because it does not have definite numerical valué 
and generally vary with the valué of x, even if C is definitised. 

The rules of integration are obtained by reversing the corre sponding 
rules of differentiation. The accuracy of this is checked by looking at the 
derivative of the integral, which must be equal to the integrand. 


Rule 1: The integral of the constant K is 
| Kdx = kx + C 

Example: J 3 dx = 3x + C 
Rule 2: The integral of 1, written as dx is 
| dx or 1 1 dx = x + C 

Rule 3: Power rule, which shows the integral of a power function 
x", where n -t -1 


f Z"dx = —-— z" +1 + C, n^-l 
J n +1 


Example 1: 


í x 2 dx = —-— 

J 2 + 1 


x 1+ ' +C = -x 2 +C 
3 


Example 2: J x 2 dx 


1 1+1 - ... 1 2 - ... 
-x +C = — x +C 

1 + 1 2 


Rule 4: Logarithmic rule: The integral of x" 1 or ]/ is, 


x 
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| x l dx = Inx + C, x > O 

The restriction x > 0 is added because only positive numbers have 
logarithms. However, for any negative numbers, it will be 

| x l dx = In\x\ + C, x 0 

Rule 5: Exponential Rule: The integral of any exponential function is 

KC 

[a^dx = —-1- C 

J teína 

2~>x 

Example: í 2 3x dx =-1- C 

J 37/72 

Rule 6: Integration of a natural exponential function 

r /** 

\£ ^dx = -1 -C, sin ce In 1 = 1 

J K 

. Q f~ 3x 

^ , Í9 r 3x dx = — — + C O' 

Example: J -3 

= -3.r 3 ; v + C 

Rule 7: Integral of a múltiple: This is the integration of a constant and an 
integrand, which is the constant (say k) times the integrated and is equal to 
the constant times the integral of the integrand. 

| Kf{x)dx = atJ f(x)dx 
Example 1: J - f(x)dx = -J f{x)dx 



3 : 1 3x 'dx = 3j x l dx = 3|n|x| + C 

Rule 8: Integral of sum or difference: The integration of sum or difference 
of two or more functions equals the sum or difference of their integráis. 

j [/(*)= g (x)\lx = J / [x)dx + 1 g ( x)dx 
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Example 1: 

| (x 3 + x +1 )dx 

f 4 X / 2 A 

= —+ C t + — + C 2 + (x + C 3 ) 

V 4 ) \ ó 


4 2 

XX 

—-1-h X + C 

4 2 


3x 3 - x + 


\)dx = 3| x 3 dx - 1 xdx + J dx 


Example 2: = 3 



—x + x + C 
2 


= —x 4 - — x 2 + x + C 
4 2 


Initial and Boundary Conditions 

The initial condition (y = yo when x = 0) or a boundary condition (y = yo 
when x = xo) uniquely determines the constant of integration. By 
permitting a unique determination of C, the initial or boundary condition 
singles out a specific curve from the family of curves. 


Example: If the boundary condition y = 11 when x = 3, the integral 


y = J 2 dx is given as 
y = 12 dx -2 x + C 


6.3 


Substituting y = 11 when x = 3, 
11 =2(3)+C 
C = 5 


Thus, y = 2x + 5. It should be noted that even though C is specified, 
12 dx remains an indefinite integral because x is unspecified. Henee, the 
integral 2x + 5 can assume an infinite number of possible valúes. 
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Integration by Substitution: The integration of a product or quotient of 
two differentiable functions of x, such as 


4 


Jl2x 2 (x 3 +2 )dx .6. 

cannot be done directly using those simple rules above. However, if the 
integrand can be expressed as a constant múltiple of another function U, 
and its derivad wedu/dx, integration by substitution is possible. Thus, 
integration by substitution is the integral of f(u)(du/dx) with respect to the 
variable x is the integral of f(u) with respect to the variable u: 

| f(u)—dx = J f{x)dx .6.5 

dx 

= J f(u)du = F(u ) + C 

Where the operation j du has been substituted for the operation j dx 

This rule is just an integral-calculus counterpart of the chain rule, which 
may be proved by chain rule itself. Given a function F(u), where u = u(x), 
the chain rules is: 

2-F(u) = ^-F(uA = f ■(»)- = m- .6.6 

dx du dx dx dx 

Given that f(u)(du/dx) is the derivative of F(u), it follows that the integral 

of the former will be; 


du 


í f(u)—dx = F(u) + C 
J dx 


6.7 


Example 1: Determine the indefinite integral of this function: 
Jl2x 2 (x 3 +2 )dx 

3 

Let u = x +2 (higher power of the independent variable in terms of 
absolute valué) 

Take the derivative of u, which is ( ^' u / c ¡ x 

Solve algebraically for dx: dx = du/ 2 

/ 

Substitute u for x 2 + 2 and du 


= 3x 


3x 


2 for dx in the original integrand: 


73 






112x^ (x 3 + 2 )dx = 112x 2 • u ■ = 1 4udu = 4jwc/w 

Where 4 is a constant múltiple of u 

= Á- U A = 2 u 2 +C = 2(x 3 + 2) 2 +C 

Example 2: find J6x 2 (x 3 + 2) 9 dx 
u = x 3 + 2; í/ze/7 = 3x 2 , so that 

16x 2 (x 3 + 2) dx = 16x 2 • u ■ = 12 udu 

3x 



= — (x 3 + 2) 10 + C 

5 V ; 


Integration by Parts 

This is applicable in a situation whereby an integrand, that is, either 
production quotient or differentiable functions of x, cannot be expressed 
as a constant múltiple of udl / dx . It is the process of reversing the product or 
quotient rule in differentiation. 

This is when an integral of v with respect to u is equal to uv less the 
integral of u with respect to v; 

| vdu = uv -1 udv .6.8 

The essence of this rule is to replace the operation J du by the 

operation J dv. The rationale behind this result is relatively simple. First, 
the product rule of differentials gives us: 
d(uv) = vdu + udv 

Integrating both sides of the equation will give 
J d (uv) = | vdu + J udv 


uv = | vdu +1 udv 
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Subtracting J udv from both sides gives previous function: 
| udv = | vdu +1 udv -1 udv 


uv — 


uv — 


j" udv = | vdu 


6.9 


Example: Intégrate by parís, this function below 

|x(x + l) 2 <ix 

Let v = x, implying dv = dx, and also let U= (X + 1) ~ so 
that du = {x + l)' l ‘ dx. Then, we can find the integral to be 
| x(x + 1 Y 2 dx = | vdu = uv -1 udv 

= ^(x + l)^x-J^(x + 1 Y 2 dx 

= — (x +1)^ x - — [x + 1)^ + C 
3 15 

Exercise: find the following 

1 . j 1 6x~ 3 dx(x ^ ()) 

2. J 4x£ x2+3 dx 

Defínite Integráis 

The theorem of fundamental calculus States that the numerical valué of the 
definite integral of a continuous function f(x) over the interval from a to b 
is given by the indefinite integral F(x) + C evaluated at the upper limit of 
integration b. minus the same indefinite integral F(x) +C evaluated at the 
lower limit of integration a. Given that C is cominon to both; the constant 
of integration is eliminated in subtraction. This is expressed 
mathematically as: 

| Íf(x)dx = F(x)\ b a = F{b)-F{a) .6.10 

| * is the limits of the integration, | b a , ]*, or [■ -] b a indicate that b and a 
are to be substituted successively for x. 


75 




Example 1: Evalúate J 10 xdx 

1 

= j“' \0xdx = 5x 2 |{= 5(4) 2 -5(1) 2 
= 80-5 = 75 

f 5 

J 3x dx AV 

Example 2: Evalúate j> 3 x 2 dx = x 3 ],=(5) 3 -(l) 3 
= 125-1 = 134 

f td x dx 

Example 3: Evalúate ^ 

\ h KB x dx=KÍ x l=(£ b -'r) 

Ja 

Properties of Definite Integráis 

Property 1: Reversing the limits of the integration changes the sign of the 
definite integral: 

í f(x)dx = - f f(x)dx 

Ja Jb 

Property 2: The valué of definite integral is zero when the two limits of 
integration are identical. 

í f (x)dx = F(a) - F(a) = 0 

Jb 

Property 3: A definite integral can be expressed as a sum of a finite 
number of definite sub-integrals as follows: 

í f{x)dx=\ f(x)dx+\ f(x)dx a<b<c 

Ja Ja Jb 

Property 4: The sum or difference of two definite integráis with identical 
limits of integration is equal to the definite integral of the sum or 
difference of the two functions. 

j f(x)dx±\ g(x)dx = f [f(x)±g(x)}jx 

Ja Ja Ja 
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Property 5: The definite integral of a constant time a function is equal to 
the constant times the definite integral of the function. 

pb pb 

rf{x)dx = ¿r f{x)dx 

Ja Ja 


Exercise: Evalúate the following 

1. f’-!v 

2 

2. f 2 (e 2x +t)dx 


Improper Integráis 

This is a situation whereby the definite integral has an infinite limit of 
integration. 

J f{x)dx and J f{x)dx . <<5v .6.11 

Are improper integráis because oois not a number and cannot be 
substituted for x in F(x). They can, however, be defined as the limits of 
other integral, such as 

p°o rb rb rb 

f(x)dx = Iim\ f{x)dx and f{x)dx=Iim\ f{x)dx .6.12 

Ja b-oo Ja J-oo a —>—co Ja 


If the limit in either case exists, then the improper integral is said to 
converge. 


Example 1: Evalúate 



dx 
Jl —yd 


dx. First we 



note 


that 


In line with the above, the desired integral is 



This improper integral does converge, and it has a valué of 1. 

The expression above is somehow cumbersome; this is the reason 
why some people prefer to write the improper integráis as: 
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f 


• dx 1 

.2 ~ 


X 


X 


= 0 + 1 


Example 2: 




í -^r- dx = lim í dx = lim 

Jl ¿7— ^oo J1 b — ^oo 


-3 


= lim 

b^>oo 


-3 (-3) 

b 1 


z' 


= Iim\ 

b—¥ ©o 


s 


— + 3 

v y 


= 3 


L’hopital’s Rule <AA 

In a situation where the limit of a function such as f(x) = g(x)/h(x) as 
x —> a cannot be evaluated, for example, when both numerator and 
denominator approach zero, giving rise to the indeterminate form 0/0 or 
when both numerator and denominator approach infinity, giving rise to the 

indeterminate form °°/°°, L' O PITAL' S rule can often be helpful. 
L’Hópital’s rule States: 


j. g(x) r <?'(+) 
lim -= Lim - 


h(x) 

Example: J¡ m 


h 1 (x) 
x-4 


.6.13 


x ^ 4 16-x 

As x —> 4, then x and 16-x 2 —>0. Therefore, differentiating 
numerator and denominator separately, we have: 

1 _ -1 
~T 


lim ———r = lim - 


- m4 16-x x ^ 4 -2x 
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LECTURE SEVEN 


Differential and Difference Equations 


Differential Equations 

These are equations that express explicitly or implicitly the relationship 
between a function y = f(t) and one or more of its derivatives or 
differentials. The following are examples of differential equations. 

* = ¿¡ + 9 . Jfe. 

dt 

.7.1 

y 1 = ]2 y 

y 11 - 2y J + 19 = 0 

These equations above are called “ordinary differential equations” 
because they have single independent variable. The integral of a 
differential equation is any equation without derivative or differential, that 
is, defined over an interval and satisfies the differential equation for all the 
valúes of the independent variable(s) in the interval. 

Example 1: to evalúate this differential equation y 11 (t) = 7 for all function 
y (t) that satisfy the equation, we simply intégrate both sides of the 
equation to find the integráis. 

yft) = jldt = lt + C l .7.2 

y(t) = j" (7f + C| )dt = 3.5 1~ + C f + C 

This is known as the general solution that shows when C is unspecified; a 
differential equation has an infinite number of possible Solutions. 


79 





However, if C is specified, the differential equation has a particular or 
definite solution, which alone of all possible Solutions is relevant. The 
differential equation order is the order of the highest derivative in the 
equation, and degree of a differential equation is the highest power to 
which the derivative of highest order is raised. 


Example 2: 

a. — = 2x + 6 
dt 


b. 


c. 


d. 


first-order, first-degree 


V dt y 

2 


-5 1 5 =0 first-order, fourth-degree 


d~y (dy 


dt 2 


+ 


\dt j 


+ = 0 second-order, first-degree 


f 1 2 A 

d y 

v dt 2 j 


+ 


v dt 3 j 


~ 75 v 


third-order, fifth-degree. 


First-Order Linear Differential Equation Formula 

The situation whereby the derivative dy/dt appears only in the first degree 
and so does the dependent variable y and whereby, no product of the form 
y (dy/dt) occurs, then the equation is called linear. The first-order linear 
differential equation will generally take the form: 

dy 


dt 


+ vy = z ■ 


.7.3 


Where v and z may be constants or functions of time, the formula for a 
general solution is: 


y(t) = £ 


-í 


vdt 


A + jz£^ dt\ .7.4 


Where A is an arbitrary constant, there are two parts to the solution: 

£^' dl | Z^' d ‘dt is called the particular integral. This particular integral y p 

equals the intertemporal equilibrium level of y(t); the complementary 
function >’r represents the deviation from the equilibrium. For y(t) to be 
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dynamically stable, y c must approach zero as t approaches infinity (that is, 
k in e kt must be negative). Solution of any differential equation can be 
checked by differentiation. 


Example: The differential equation general solution of this function dy/dt 


+ 4y = 12 is given as follows: since v = 4 and z = 12, substitute in 
equation (7.4) 

Y(t) = Ó 4dt j^A + J 12.^ 4 '" dt j.7.5 

From the integration rules, J4 dt = 4t + C. In equation (7.5) C is ignored 
and subsumed under A. Thus, 

y(t) = £~ 4 '(á + jl2£ 4t dt) ...7.6 

Integrating the remaining integral gives j 1 21 4l dt = 37 4 ' + C . Ignoring C 
again and substitute in equation (7.6): 

y(t) = V* (A + ) = AC At +3.7?...7.7 

Since 


1~ 4, £ At = £° = 1. As t —y °°, y p = Ar M —> 0 and y(t) approachs y p = 3, 

the intertemporal equilibrium level. y(t) is dynamically stable. All these 
illustrations above are what are known as the non-homogenous case of 
differential equations. 


The Homogenous Case 

Given the function — + u(t)y = w(t) .7.8 

dt 


If u and w are constant functions and if w happens to be identically zero, 
then equation (7.8) will become: 


dy 

— + ay = 0. 
dt - 


.7.9 


Where a is constant, this differential equation is homogenous on the 
account of the zero constant term. The defining features of a homogenous 
equation is that when all the variables (here dy/dt and y) are multiplied by 
a given constant, the equation remain valid. This feature holds if the 
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constant term is zero, but will be lost if the constant term is not zero. An 
alternative to equation 2.2 is given as: 


1 dy 

- - = -a 

y dt 


.7.10 


The solution to equation (7.9) or (7.10) is 

y(t) = Ai~ at general solution.7.11 

Or y(t) = y(Q¡)C at definite solution.7.12 


General solution occurs as a result of the arbitrary constant A, but when 
any particular valué is substituted for A, the solution becomes a particular 
solution of equation (7.9). There is an infinite number of particular 
Solutions, one for each possible valué of A, including the valué y (0), 
which is the only valué that can make the solution satisfy the initial 
condition. Equation (7.12) will be referred to as the definite solution of 
the differential equation (7.9) or (7.10). Since it represents the result of 
definitising the arbitrary constant. 

Exercise 

2 

1. find the general solution of dy/dt + Sry 
t 2 

2. Calcúlate the definite solution when y(0) 

Second-Order Differential Equations 

These equations require sepárate Solutions for the complementary function 
y c and the particular integral y p . The general solution is the sum of the two: 
y(t) = y- r + y p . Given the second-order linear differential equation. 

T(t) + b iy¡ (t) + b 2 y(t) = a .7.13 

Where b¡, b 2 and a are constants, the particular integral will be, 


2 2 

= t, where v = 3t and z = 

= 1 in question (1). 





b 2 * 0 


7.14 



¿ 2 = 0 , ^ 0 


.7.15 
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y p =-t 2 \=b 2 = O.7.16 

The complementary function is 
Where 

= yi + y 2 . 7 - 17 




And r,, r 2 =■ 


b \ 


. .O.....7.20 


Here, A; and A 2 are arbitrary constants, and Z> 2 -t- 4b 2 

r¡ and r 2 are referred to as characteristic roots, and equation (7.20) is the 
solution to the characteristic or auxiliary equation: r + b¡r + = 0. 


Example 1: Find the particular integral for each of the following 
equations: 


i. y 11 <t)-5y(t) + 4y(t) = 2 


2 1 


using equation 30, we have y p - — = —.7.21 


ii. y n (t) + 3y ! (t) = 12 

c-> 12 

using equation 31, y p = — t = 41 .7.22 

iii. y u (t) = 7(5 

<Y\~ 1 6 T T 

using equation 32, y = — t 2 = 8r.7.23 


Example 2: The complementary functions for equation 1 and 2 in example 
1 above are calculated as: 

„ , +5±V(-5) 2 -4(4) 5 ± 3 

For equation 1: rjr 2 =--—--= = 1,4 

Substituting in equations (3.1) and (3.2) and finally in (3.0) 
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y c =A/+A 2 l At 


7.24 


c „ -3±V(3) 2 -4(0) -3 + 3 

For equation 2 r t r 2 =---= —-— = 


Thus, y c = A/° + A 2 r 3í = A, + A 2 f 3í .7.26 

= 0±V0-4(0) = 

For equation 3: ^ 2 .2.27 


y c =A/ 0 +V 0 =A 1+ A 2 

Example 3: The general solution of a differential equation is composed of 
the complementary function and the particular integral, that is, y(t) = y c + 
y p . From example 1, the general Solutions to equation (1) and (2) are: 


For equation 1, x 7\ 

y(t) = A/ 1 + A 2 t' + l / 2 ..7.28 

For equation 2, 

y(t) = A, + A 2 r 3t +4 1 .7.29 


Example 4: The definite solution for equation (4.0) is calculated as; let 
5 K and Ny (0) — 11, then, 

assume y(0) = 5 .7.30 

y(t) = A l t+A 2 t'+y 2 

Thus, y\t) = AA\ + 4AJ 4t .7.31 


Evaluating equation (4.2) and (4.3) at t = 0, and setting y(0) = 
5 x / 2 and y 1 (0) = 1 from the initial conditions, 

d$j=A 1 £°+A 2 e 4m +y 2 =5y 2 

= A 1 + A 2 =5 .‘' , " 

^‘(O) = A/° +4A 2 ( 4<0) =11 7 „ 

= ,\. l 4A. = I 1 . ' " 

Solving equation (4.4) and (4.5) simultaneously, A¡ = 3 and A 2 = 2, 
substituting in equation (4.0) 

y(t) = 3f +2£ 4 ’ +y .7.34 
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Difference Equations 

The difference equations are equations that show the relationship between 
a dependent variable and lagged independent variables, which changes at 
discrete interval of time. For example, C t = /'(y, j), where I, Y and T 
are measured at the end of each year. The difference equation order is 
determined by the greatest number of periods lagged. The difference 
equation first-order shows a time lag of one period; a second-order, two 
periods. The change in y as t changes from t to tt, is called the first 
difference of y, which is written as: 

~ - Ay = y ¡ - y .,47.35 

At 

The pattern of change is now represented by the difference quotient 
Ay/At, which is the discrete-time counterpart of the derivative dy/át. The 
solution of a difference equation defines y for every valué of t and does 
not contain a difference expression. 

Example 1: 


+ 

II 

1 

order 1 

Qs t -1 =a + bp t _ 2 

order 2 

Y f+3 -9y í+2 + 2 y f+ i + 6y ? = 8 

order 3 

Ay, = 5 y, 

order 1 


Example 2: If the initial valué of y is yo, the difference equation 

y (A = by t .7.36 

Solution by successively substituting t = 0,1,2,3, etc will be 
?= 0 y, = by 0 

t = 1 =by x =b(by 0 ) =b 2 y 0 

t = 2 y 3 =by 2 =b(b 2 y 0 ) = b i y 0 

t = 3 y 4 = by 3 = b(b 3 y 0 ) = b*y 0 
Therefore, for any period, 
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y, =b'y 0 

This method is called the iterative method. Since yo is a constant, notice 
the crucial role b plays in determining valúes for y as t changes. 


First-Order Linear Difference Equations General Formula 

A linear first-order difference equation is an equation where all the 
variables are raised to the first power, without any cross products. For 
example, 

y, =by t -i+a .:4\ 7 - 37 

Where b and a are constants, the general formula for a definite solution is 

( ~ \ 


y t 


Jo 


V 


a 

I ~b 


\b' + ■ 


J 


1 -b 


when b -t 1.7.38 


y t = y 0 + at 


when b = 1 


.7.39 


In a case where no initial condition is given, an arbitrary constant A is 
used for yo - a/(l-b) in equation (7.38) and for yo in equation (7.39). This 
is called the general solution. 

Example: Solve this difference equation. 

y t = -7 y,_ x +16 and y 0 = 5. In the equation, b = -7 and a = 16 
b^Y, then 


Since 


y¡ 


16 7)' + 16 


1 + 7 


.7.40 


1 + 7 


= 3(-7 y + 2 

To cheek the answer, substitute t = 0 and t = 1 in equation (7.37) 

Jo =3(-7)° +2-5, sin ce (-7)° = 1 
Ji = 3(-7) : + 2 = -19 

Substituting yi = -19 for y t and yo = 5 for y t _i in the original equation, 
-19 =-7(5)+ 16 = -35 + 16. 
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The Dynamic Stability Conditions 

In general form, equation (7.38) can be expressed as: 
y t = Ab' + C 

In equation (7.40) A = yo - a/(l - b) and C = a/(l - b). Ab' in equation 
(54) is called the complementary function and c is the particular solution. 
The intertemporal equilibrium level of y is expressed by the particular 
solution while the complementary function represents the deviations from 
that equilibrium. Equation (7.40) will be dynamically stable only if the 
complementary function Ab' —> 0, as t —> °° . All these will depend on 
the base b. Assuming A = 1 and C = 0 for the moments, the exponential 
expression b 1 will generate seven different time paths depending on the 
valué of b. 
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Example: Time path of b*. 
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From these diagrams, if \b\ > 1, the time path will explode and move farther 
and farther away from equilibrium, if|¿>| < 1, the time path will be damped 

and move toward equilibrium. If b < 0, the time path will oscillate 
between positive and negative valúes, if b > 0, the time path will be non- 
oscillating. If A ^ 1, the valué of the multiplicative constant will scale up 
or down the magnitude of b 1 . but will not change the basic pattern of 
movement. If A = -1, a mirror image of the time path of b 1 with respect to 
the horizontal axis will be produced. If C ^ 0, the vertical intercept of the 
graph is affected, and the graph shifts up or down accordingly. 

Exercise 

(a) solve the difference equation given below; (b) check your answer, 
using t = 0 and t = 1; and (c) comment on the nature of the time path. 

y, = 6y 

5y f +2y f _i-140 = 0 y 0 =3 ° 


Second-Order Difference Equations 

The second-order difference equations general solution composed of both 
complementary function and particular solution: y(t) = y c + y p . Given the 
second-order linear difference equation. 

^ + b J t -i 1 b ?y, : = «.7.41 

Where bi, b 2 , and a are constants, the particular solution is: 


1 + Z?i + b 2 


b \ + b 2 ^ ~ 1 


7.42 


y = t b. +b 2 =-l b,*- 2.7.43 

y 2 + b x 


y p = ^r b l+ b 2 =- 1 b,=- 2.7.44 

and the complementary function is: 

A: = V/ + A 2 r 2 .7.45 
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Where A¡ and A 2 are arbitrary constants and the characteristic roots r¡ and 
r 2 are found using 

b\ ± V¿>i 2 ~ 4¿> 2 
2 

Given that bf ^ 4 b 2 


Example 1: Find the particular solution for the following equations 

1- 3",-10y^ +16y f _ 2 =14 ^ 

2- y, -6y,.i +5y f _ 2 =12 

3 y, - 2y t _ l + y,_ 2 = 8 

For question (1), using equation (7.42) 

14 

y =-= 2 

p 1-10 + 16 

12 -t = -3t 


2-6 


89 . 2 

>’- 2 ‘ 


Example 2: The complementary function of equation (1) and (2) in 
example 1 using equation (33) and then substitute in equation (7.45) is: 

101-^100-4(16) 1016 

r.r 2 =-=-= 2,8 

,2 2 

thus, y c = Aj(2)' + A 2 (8)' 

for equation(2) 

6± ^36-4(5) 614 

= -=-= 1,5 

“ 2 2 

thus y c = A,(l)'+ A,(5)' = A 1 +A 2 (5) f 


Characteristic Roots 

The characteristic equation have three different types of roots, they are: 
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1. Distinct real roots: Given thatZ?/ > 4Z?,, the square root in equation 
(33) will be a real number, and ri and n will be distinct real 
numbers. 


2. Real roots that are repeated: If Z?, 2 = 4 b 2 , the square root of 
equation (33) will vanish, and ri and r 2 will equal the same real 
number. The formula for y c in repeated real roots will be: 


y c = A v r + A 2 +r 


.7.46 


y c - A r? + A + r ' .7.47 

3. Complex roots: Equation (33) will contain the square root of a 
negative number if bf < 4b ,, and this is called imaginary number. 
In this case, ri and tt are complex numbers. A complex number 
contains a real part and an imaginary part; for example, (12 +i) 
where i = 


Example: the complementary function for equation 3 in example 1 above, 
where y n (t) = 16 is 

0±V0 + 4(0) 

r,r 2 =-= 0 

2 

Based on equation (60), Tu Xz = 0, which is a case of repeated real roots, 
y c — A/° + A 2 + £° = Aj + A 2 t. In question (3) of example, 

37 - 2y ; _, + y,^ = 8 

Solving for the complementary function, from equation (33), 

"2 ± a/4-4(1) 2±0 , 

r.r^ = -=-= 1 

1 <^r 2 2 

However, based on equation (61), since ri = r 2 = 1, 
y c =A(l)'+A 2 í(l)'=A 1+ A 2 t 
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LECTURE EIGHT 


Optimal Control Theory 


Introduction 

The optimal control theory is an advanced methodology to the classical 
approach to dynamic optimization called calculus of variations. This 
theory can handle any problem the calculus of variations was designed for, 
and further manage some problems that the calculus of variation cannot, 
such as córner solution, and constraints on the derivatives of the functions 
sought. 

The basic aim of this theory is to find the optimal time path for a 
control variable, which can be denoted as y. In the calculus of variations, 
the variable with which we sought an optimal time path is known as ‘a 
State variable’, which we shall designate as x. The State variables often 
have equations of motion set equal to x. Thus, the goal of the theory of 
optimal control is to select a stream of valúes over time for the control 
variable that will optimize the functional subject to the constraint set on 
the State variable. 

The optimal control theory problems that involve continuous time, a 
finite time horizon, and fixed endpoint are generally written as: 


Maximize 

/=[ f[x{t),y(t),t]it .8.1 

JO 

Subject to x = g[x{t),y{t),t] . 8 . 2 


x(0) = x 0 x{T) = x T 

In equation (8.1), J is the valué of the function to be optimized, y(t) is the 
control variable, because its valué is selected to optimize. J; x(t) is the 
State variable that changes over time according to the differential equation 
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set equal to xin the constraint, and whose valué is indirectly determined 
by the control variable in the constraint; and t is the time. The solution to 
an optimal control problem demarcates the optimal dynamic time path for 
the control variable y(t). 

Example: Lets consider an economy that produces output Y using capital 
K and a fixed amount of labour L, the production function of the economy 
is given as: 

Y=Y(K,L) 

Output can be used either for consumption C or investment I. if we ignore 
the problem of depreciation, then 

_ dK 
dt 

This means that investment is the change in capital stock overtime. Thus, 
we can also write investment as: 

dK 

I = Y-C = Y(K,L)-C = - 

dt 

which gives us a first-order differential equation in the variable K. 

If the objective is to maximize some form of social utility over a fixed 
planning period, then the problem becomes 

Maximise í U ( C)dt 

Jo 

Subject to < ^~= Y(K,L) — C 

dt 

And K(0) = K 0 K(T) = K T 

Where Ko and Kj are the initial valué and terminal valué of K, 
respectively. Note that in the above example, the terminal State is a fixed 
valué and not left free. Here C serves as the control variable and K is the 
State variable. The problem is to choose the optimal control path C(t) such 
that its impact on output Y and capital K, and the repercussions therefrom 
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upon C itself, will together maximize the aggregate utility over the 
planning period. 


The Máximum Principie 

This principie is attributed to L.S. Pontryagin and his associates, and is 
often referred to as Pontryagin’s máximum principie. This is the first- 
order necessary condition of the optimal control theory. The statement of 
the máximum principie involves an approach that is akin to the 
Lagrangian function and the Lagrangian multiplier variable. For the 
optimal control problems, these are known as Hamiltonian function and 
cosíate variable. The dynamic optimization of a function subject to a 
constraint on the State variable in optimal control involves a Hamiltonian 
function H similar to the Lagrangian function in concave programming. 
With respect to equation (62), the Hamiltonian is defined as: 

H[x(t), y(t), A(t), t] = f[x(t), y(t),t]+ A{t)g[x{t), y(t), t\ .8.2 


Where A(í) is called the cosíate variable. Similar to the LaGrangian 
multiplier, the cosíate variable A(t) estimate the marginal valué or shadow 
price of the associated State variable x(t). In equation (8.2), the formation 
of the Hamiltonian is easy and this is by taking the integrand under the 
integral sign and adds to it the product of the costate variable A(t) times 
the constraint. 


Let’s assume that the Hamiltonian is differentiable in y and strictly 
concave so that there is an interior solution and not an endpoint solution, 
the necessary conditions for maximization are 

. SH n 


2 (a) 2 3 * S A = A = ^- 

& Sx 

3 (a) x(0) = x 0 


& * SH 

ib) — = x =- 

& &A 

( b ) x(T) = x T 


The first two conditions are known as the máximum principie and the 
third is called the boundary condition. The two equations of motion in the 

second condition are generally referred to as the Hamiltonian system or 
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the Canonical system. To maximize, the objective function can simply be 
multiplied by -1, as in the concave programming. And if the solution does 
not involve an end point, the, SH / Sy need not equal to zero in the first 
condition, however, H must be maximized with respect to y. We shall 
generally assume interior Solutions. 


Example: Using the conditions in máximum principie, 

Maximize [ (4.x - 5 y 2 )dt 

jo 

Subject to * = 

x(0) = 2 x(3) = 117.2 

A. Set up the Hamiltonian 

H = 4x-5y 2 + >1(8 y) 

B. Assume an interior solution by applying the máximum principie 

-=0 

SH 

1. — = -10y + 8/l = 0,.:.:..8.3 

oy 

y = 0.8/1 

Á = ~ — 

2(a) Sx .8.4 

^ í = -4 
' _ SH 
(b) * SA 

x = -8 y 

From equation (8.3), y = 0.8/1, then, 

x = S(0.M) = 6AA .8.5 
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Since we have employed the máximum principie, we are now left with 
two differential equations, which we shall solve for their State variables 
x(t) and the costate variable A(t). 

Integrating equation (8.4) will give us the State variable 

Á{t) = \\dt = \-Adt = -At + C, .8.6 

Substitute equation (8.6) and (8.5) 

x = 6.4(-4í + C,) = -25.6 1 + 6.4Q 
Integrating, 

x(t) = í(- 25.6 1 + 6.4C, )dt 

J ..8.7 

x(t) = -12.8í 2 + 6.4C,t + C 2 

C. To solve for the constants of integration, we use the boundary 
conditions where x(0) = 2, x(3) = 117.2 and apply to equation (8.7) 

x( 0) = -12.8(0) 2 +6.4C 1 (0) + C 1 =2 C 2 = 2 
x(3) = -12.8(3) 2 + 6.4C, (3) + 2 = 117.2 C 2 = 12 
Thus, substituting Cj = 12 and C 2 = 2 in equation (68) and (67), we have; 

x(t) = 12.8t 2 +6.4(12)t+2 .8.8 

x(t) = 12.8t 2 + 76.8t + 2 State variable 
/l(0 = -4í + 12 costate variable 

D. The solution for the control variable y(t) can be derived in two ways: 

1. In equation (8.3), y(t) = 0.8/1, so 

y(t) = 0.8(-4t+12) = -3.2t + 9.6 control variable 

2. Or take the derivative of equation (8.8) 

x = -25.6 + 76.8 
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Then substitute for xin the equation of motion in the constraint, 
x = 8 v 

- 25.6 1 + 76.8 = 8 y 

y(t ) = -3.2 1 + 9.6 Control variable 

Evalúate at the endpoints, we have, 
y(0) = -3.2(0) + 9.6 = 9.6 
y(3) =-3.2(3)+ 9.6 = 0 

Therefore, the optimal path of the control variable is linear starting at 
(0,9.6) and ending at (3,0), with a slope of -3.2. 

Sufficient Conditions 

Assuming that the necessary condition which is the máximum principie, 
for maximization in optimal control is satisfied, ten the sufficient 
conditions will be fulfilled if: 

1. Both the objective functional f[x(t), y(t), t] and the constraint 
g[x(t), y(t),t] are differentiable and jointly concave in x and y, and; 

2. A(t) > 0, if the constraint is nonlinear in x or y. If the constraint in 
linear, the A may assume any sign. 

Linear functions are always both concave and convex, but neither 
strictly concave ñor strictly convex. However, an easy test for joint 
concavity is the simple discriminant test for non-linear functions. The 
discriminant of the second-order derivatives of a function, 

D = fax fxy 

fyx fyy 

Will be strictly concave if the discriminant is negative definite, 

\D t \ = fxx<0 and |D 2 | = |d|> 0 
And simply concave if the discriminant is negative semi-definite, 

|Dj | = fxx < 0 and | D 2 | = |d| > 0 

A negative definite discriminant indicates a global máximum and is, 
therefore, always sufficient for a máximum. An indicative of a local 
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máximum is the negative semi definite, which is sufficient for a 
máximum, if the test is conducted for every possible ordering of the 
variables with similar results. 


Example: The sufficient conditions for the problem in the example under 
the máximum principie are given below. Let us start with the objective 
functional that is non-linear, we take the second derivatives and apply the 
discriminant test. 


D = 


fxx 

Jyx 


fxy 

fyy 


o o 

o -10 


where D, = 0 and D 2 = D = 0 


Thus, since both discriminant tests are negative semi-definite, then 
objective functional f is jointly concave in x and y. Given that the 
constraint is linear, it is also jointly concave and does not need testing. 
Therefore, we conclude that the functional is indeed maximized. 
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